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Abstract 

We consider a model of half-filled bilayer graphene, in which the three domi¬ 
nant Slonczewski-Weiss-McClure hopping parameters are retained, in the presence of 
short range interactions. Under a smallness assumption on the interaction strength 
U as well as on the inter-layer hopping e, we construct the ground state in the ther¬ 
modynamic limit, and prove its analyticity in U, uniformly in e. The interacting 
Fermi surface is degenerate, and consists of eight Fermi points, two of which are 
protected by symmetries, while the locations of the other six are renormalized by 
the interaction, and the effective dispersion relation at the Fermi points is coni¬ 
cal. The construction reveals the presence of different energy regimes, where the 
effective behavior of correlation functions changes qualitatively. The analysis of the 
crossover between regimes plays an important role in the proof of analyticity and in 
the uniform control of the radius of convergence. The proof is based on a rigorous 
implementation of fermionic renormalization group methods, including determinant 
estimates for the renormalized expansion. 
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1 Introduction 


Graphene, a one-atom thick layer of graphite, has captivated a large part of 
the scientific community for the past decade. With good reason: as was shown by 
A. Geim’s team, graphene is a stable two-dimensional crystal with very peculiar elec¬ 
tronic properties [NGe04]. The mere fact that a two-dimensional crystal can be synthe¬ 
sized, and manipulated, at room temperature without working inside a vacuum [GelO] 
is, in and of itself, quite surprising. But the most interesting features of graphene 
lay within its electronic properties. Indeed, electrons in graphene were found to have 
an extremely high mobility [NGe04], which could make it a good candidate to replace 
silicon in microelectronics; and they were later found to behave like massless Dirac 
Fermions [NGe05, ZTe05], which is of great interest for the study of fundamental Quan¬ 
tum Electro-Dynamics. These are but a few of the intriguing features [GN07] that have 
prompted a lively response from the scientific community. 

These peculiar electronic properties stem from the particular energy structure of 
graphene. It consists of two energy bands, that meet at exactly two points, called the 
Fermi points [Wa47]. Graphene is thus classified as a semi-metal: it is not a semi¬ 
conductor because there is no gap between its energy bands, nor is it a metal either 
since the bands do not overlap, so that the density of charge carriers vanishes at the 
Fermi points. Furthermore, the bands around the Fermi points are approximately con¬ 
ical [Wa47], which explains the masslessness of the electrons in graphene, and in turn 
their high mobility. 

Graphene is also interesting for the mathematical physics community; its free 
energy and correlation functions, in particular its conductivity, can be computed non- 
perturbatively using constructive Renormalization Group (RG) techniques [GMIO, GMPll, 
GMPI2], at least if it is at half-filling, the interaction is short-range and its strength is 
small enough. This is made possible, again, by the special energy structure of graphene. 
Indeed, since the propagator (in the quantum field theory formalism) diverges at the 
Fermi points, the fact that there are only two such singularities in graphene instead 
of a whole line of them (which is what one usually finds in two-dimensional theories), 
greatly simplifies the RG analysis. Furthermore, the fact that the bands are approxi¬ 
mately conical around the Fermi points, implies that a short-range interaction between 
electrons is irrelevant in the RG sense, which means that one need only worry about the 
renormalization of the propagator, which can be controlled. 

Using these facts, the formalism developed in [BG90] has been applied in [GMIO, 
GMPI2] to express the free energy and correlation functions as convergent series. 

Let us mention that the case of Goulomb interactions is more difficult, in that 
the effective interaction is marginal in an RG sense. In this case, the theory has been 
constructed at all orders in renormalized perturbation theory [GMPIO, GMPIlb], but 
a non-perturbative construction is still lacking. 

In the present work, we shall extend the results of [GMIO] by performing an RG 
analysis of half-filled hilayer graphene with short-range interactions. Bilayer graphene 
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consists of two layers of graphene in so-called Bernal or AB stacking (see below). Before 
the works of A. Geim et al. [NGe04], graphene was mostly studied in order to understand 
the properties of graphite, so it was natural to investigate the properties of multiple 
layers of graphene, starting with the bilayer [Wa47, SW58, Mc57]. A common model 
for hopping electrons on graphene bilayers is the so-called Slonczewski-Weiss-McClure 
model, which is usually studied by retaining only certain hopping terms, depending on 
the energy regime one is interested in; including more hopping terms corresponds to 
probing the system at lower energies. The fine structure of the Fermi surface and the 
behavior of the dispersion relation around it depends on which hoppings are considered 
or, equivalently, on the range of energies under inspection. 

In a first approximation, the energy structure of bilayer graphene is similar to 
that of the monolayer: there are only two Fermi points, and the dispersion relation is 
approximately conical around them. This picture is valid for energy scales larger than 
the transverse hopping between the two layers, referred to in the following as the first 
regime. At lower energies, the effective dispersion relation around the two Fermi points 
appears to be approximately parabolic, instead of conical. This implies that the effective 
mass of the electrons in bilayer graphene does not vanish, unlike those in the monolayer, 
which has been observed experimentally [NMe06]. 

From an RG point of view, the parabolicity implies that the electron interactions 
are marginal in bilayer graphene, thus making the RG analysis non-trivial. The flow 
of the effective couplings has been studied by O. Vafek [ValO, VYIO], who has found 
that it diverges logarithmically, and has identified the most divergent channels, thus 
singling out which of the possible quantum instabilities are dominant (see also [TV12]). 
However, as was mentioned earlier, the assumption of parabolic dispersion relation is 
only an approximation, valid in a range of energies between the scale of the transverse 
hopping and a second threshold, proportional to the cube of the transverse hopping 
(asymptotically, as this hopping goes to zero). This range will be called the second 
regime. 

By studying the smaller energies in more detail, one finds [MF06] that around 
each of the Fermi points, there are three extra Fermi points, forming a tiny equilat¬ 
eral triangle around the original ones. This is referred to in the literature as trigonal 
warping. Furthermore, around each of the now eight Fermi points, the energy bands 
are approximately conical. This means that, from an RG perspective, the logarithmic 
divergence studied in [ValO] is cut off at the energy scale where the conical nature of the 
eight Fermi points becomes observable (i.e. at the end of the second regime). At lower 
energies the electron interaction is irrelevant in the RG sense, which implies that the flow 
of the effective interactions remains bounded at low energies. Therefore, the analysis of 
[ValO] is meaningful only if the flow of the effective constants has grown significantly in 
the second regime. 

However, our analysis shows that the flow of the effective couplings in this regime 
does not grow at all, due to their smallness after integration over the first regime, which 
we quantify in terms both of the bare couplings and of the transverse hopping. This puts 
into question the physical relevance of the “instabilities” coming from the logarithmic 
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divergence in the second regime, at least in the case we are treating, namely small 
interaction strength and small interlayer hopping. 

The transition from a normal phase to one with broken symmetry as the interac¬ 
tion strength is increased from small to intermediate values was studied in [CTV12] at 
second order in perturbation theory. Therein, it was found that while at small bare cou¬ 
plings the infrared flow is convergent, at larger couplings it tends to increase, indicating 
a transition towards an electronic nematic state. 

Let us also mention that the third regime is not believed to give an adequate 
description of the system at arbitrarily small energies: at energies smaller than a third 
threshold (proportional to the fourth power of the transverse hopping) one finds [PP06] 
that the six extra Fermi points around the two original ones, are actually microscopic 
ellipses. The analysis of the thermodynamic properties of the system in this regime (to 
be called the fourth regime) requires new ideas and techniques, due to the extended 
nature of the singularity, and goes beyond the scope of this paper. It may be possible 
to adapt the ideas of [BGM06] to this regime, and we hope to come back to this issue 
in a future publication. 

To summarize, at weak coupling and small transverse hopping, we can distinguish 
four energy regimes: in the first, the system behaves like two uncoupled monolayers, 
in the second, the energy bands are approximately parabolic, in the third, the trigonal 
warping is taken into account and the bands are approximately conical, and in the fourth, 
six of the Fermi points become small curves. We shall treat the first, second and third 
regimes, which corresponds to retaining only the three dominant Slonczewski-Weiss- 
McClure hopping parameters. Informally, we will prove that the interacting half-filled 
system is analytically close to the non-interacting one in these regimes, and that the 
effect of the interaction is merely to renormalize the hopping parameters. The proof 
depends on a sharp multiscale control of the crossover regions separating one regime 
from the next. 

We will now give a quick description of the model, and a precise statement of the 
main result of the present work, followed by a brief outline of its proof. 


1.1 Definition of the model 

We shall consider a crystal of bilayer graphene, which is made of two honeycomb 
lattices in Bernal or AB stacking, as shown in figure 1.1. We can identify four inequiva¬ 
lent types of sites in the lattice, which we denote by a, b, d and b (we choose this peculiar 
order for practical reasons which will become apparent in the following). 

We consider a Hamiltonian of the form 

n = nQ + ni (i.i) 

where the free Hamiltonian Ho plays the role of a kinetic energy for the electrons, and 
the interaction Hamiltonian Hi describes the interaction between electrons. 
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Figure 1.1; 9 and ■ represent atoms of type a and b on the lower layer and Q and 
I I represent atoms of type a and b on the upper layer. Full lines join nearest neighbors 
within the lower layer and dashed lines join nearest neighbors within the upper layer. 


Tio is given by a tight-binding approximation, which models the movement of 
electrons in terms of hoppings from one atom to the next. There are four inequivalent 
types of hoppings which we shall consider here, each of which will be associated a different 
hopping strength 7 *. Namely, the hoppings between neighbors of type a and b, as well as 
a and b will be associated a hopping strength 70 ; a and b a strength 71 ; a and b a strength 
73 ; a and a, and b and b a strength 74 (see figure 1.2). We can thus express Hq in second 
quantized form in momentum space at zero chemical potential as [Wa47, SW58, Mc57] 

= ( 1 . 2 ) 

' fceA 


M := 


and Hq^Iv) := — 


A 71 74^2(/c) 

71 A 7o^2(/c) 

7417* (/c) 7 o17*(/c) 0 

7o 17(A:) 74l7(A:) {k)e~^^^^ 


70^* (k) 
7417* (A:) 
73l7(/c)e^*^^ 


in which dk, bk, dt and bk are annihilation operators associated to atoms of type a, b, d 
and b, k = {kx,ky), A is the first Brillouin zone, and 17(/c) := 1 + cos 

These objects will be properly defined in section 2.1. The A parameter in Hq models a 
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shift in the chemical potential around atoms of type a and b [SW58, Mc57]. We choose 
the energy unit in such a way that 70 = 1. The hopping strengths have been measured 
experimentally in graphite [DD02, TDD77, MMD79, DDe79] and in bilayer graphene 
samples [ZLe08, MNe07]; their values are given in the following table: 


(1.4) 


We notice that the relative order of magnitude of 73 and 74 is quite different in graphite 
and in bilayer graphene. In the latter, 71 is somewhat small, and 73 and 74 are of the 
same order, whereas A is of the order of 7 ^. We will take advantage of the smallness of 
the hopping strengths and treat 71 =: e as a small parameter: we fix 

II = 1 , 71 = 0 . 33 , II = 0.40, ^ = 0.58 (1.5) 

6 6 6 6 ^ 

and assume that e is as small as needed. 

Remark: The symbols used for the hopping parameters are standard. The reason why 
72 was omitted is that it refers to next-to-nearest layer hopping in graphite. In addition, 
for simplicity, we have neglected the intra-layer next-to-nearest neighbor hopping 7 q ~ 
O.I 71 , which is known to play an analogous role to 74 and A [ZLe08]. 

The interactions between electrons will be taken to be of extended Hubbard form, 

i.e. 

1 
2 

where Ux ■= extc^x in which ax is one of the annihilation operators bx, dx or bx', the 
sum over (x, y) runs over all pairs of atoms in the lattice; u is a short range interaction 
potential (exponentially decaying); U is the interaction strength which we will assume 
to be small. 





We then define the Gibbs average as 

(■) := Ixr 

where 

Z := Tr =: 

The physical quantities we will study here are the free energy /, and the two-point 
Schwinger function defined as the 4x4 matrix 

.S2(xi-X 2 ) := ((T(a' ai )\) (1.7) 

\ \ / / {a' ,<y)£{a,b,a,b}^ 
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Figure 1.2: The different types of hopping. From top-left to bottom-right: a -H- 6, a -H- 6, 
a b, b a, a a, b ^ b. Atoms of type a and a are represented by spheres and 
those of type b and b by cubes; the interaction is represented by solid red (color online) 
cylinders; the interacting atoms are displayed either in purple or in blue. 
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where xi := and X 2 ;= (^ 2 )^: 2 ) includes an extra imaginary time component, 

ti ,2 £ [0,/3), which is introduced in order to compute Z and Gibbs averages, 

C(t,x '■= for a G {a, b, a, b} 


and T is the Fermionic time ordering operator. 


T(a' 




if > ^2 

-«i2,X2«tl,Xl if ^ ^2 


( 1 . 8 ) 


We denote the Fourier transform of S 2 (x) (or rather of its anti-periodic extension in 
imaginary time for t’s beyond [0, j3)) by S 2 (k) where k := (/cq. A;), and /cq £ ^(^ + 5 )- 


1.2 Non-interacting system 

In order to state our main results on the interacting two-point Schwinger function, 
it is useful to first review the scaling properties of the non-interacting one, 

=-{ikot + Ho{k))-\ 

including a discussion of the structure of its singularities in momentum space. 

1 - Non-interacting Fermi snrface. If HQ{k) is not invertible, then s®(0,A:) 
is divergent. The set of quasi-momenta Tq := {k, det HQ{k) = 0} is called the non¬ 
interacting Fermi surface at zero chemical potential, which has the following structure; 
it contains two isolated points located at 

Ppfi •= ^ ^ {“f’+f} (f-9) 

around each of which there are three very small curves that are approximately elliptic 
(see figure 1.3). The whole singular set J-q is contained within two small circles (of 
radius 0 (e^)), so that on scales larger than e^, Tq can be approximated by just two 
points, figure 1.3. As we zoom in, looking at smaller scales, we realize that 

each small circle contains four Fermi points: the central one, and three secondary points 
around it, called {Ppj, j £ {1)2,3}}. A finer zoom around the secondary points reveals 
that they are actually curves of size e^. 

2 - Non-interacting Schwinger function. We will now make the statements 
about approximating the Fermi surface more precise, and discuss the behavior of 
around its singularities. We will identify four regimes in which the Schwinger function 
behaves differently. 

2-1 - First regime. One can show that, if p^Q ;= (0,p^g), and 

||(A;o,A;',fe;)||i := ^k^ + {k'^f + {k'yf 
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Figure 1.3: Schematic representation of the Fermi points. Each dotted square represents 
a zoom into the finer structure of the Fermi points. The secondary Fermi points are 
labeled as indicated in the figure. In order not to clutter the drawing, only one of the 
zooms around the secondary Fermi point was represented. 












Figure 1.4: Eigenvalues of Ho{k). The sub-figures b,c,d are finer and finer zooms around 
one of the Fermi points. 


then 

4°^(Pf,o+ 1^0 = ('2i^(Pf,o+ 1 ^ 0 ) ( 1 + 0(||k'||i,e||k'||f^)) (1.10) 

in which £,iA is a matrix, independent of 71 , 73 , 74 and A, whose eigenvalues vanish 
linearly around p^Q (see figure 1.4b). We thus identify a first regime: 

e <C ||k'||i <C 1 

in which the error term in (1.10) is small. In this first regime, 71 , 73 , 74 and A are 
negligible, and the Fermi surface is approximated by {p^'q}, around which the Schwinger 
function diverges linearly. 

2-2 - Second regime. Now, if 

\\{ko,k'^,ky)\\u := ^k"^ + 

then 

4°^(Pf,o + ^0 = ('2ii^(Pf,o + k')) ( 1 + 0(e“^||k'||ii,e^/2||k'||"//^)) (1.11) 

in which S^uA is a matrix, independent of 73 , 74 and A. Two of its eigenvalues vanish 
quadratically around p^Q (see figure 1.4c) and two are bounded away from 0. The latter 
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correspond to massive modes, whereas the former to massless modes. We thus identify 
a second regime: 

< llk'llii < e 

in which 73 , 74 and A are negligible, and the Fermi surface is approximated by {p^q}, 
around which the Schwinger function diverges quadratically. 

2-3 - Third regime. If := j = 0,1, 2, 3, and 

\\{ko,k'j^^,kj^y)\\ui := \Jkl + 

then 

+ k')) (l + 0(e“3||k'.||iii, e^||k'.||"j^)) (1.12) 

in which S^m^jA is a matrix, independent of 74 and A, two of whose eigenvalues vanish 
linearly around ;= {0,p^j) (see figure 1.4d) and two are bounded away from 0. We 
thus identify a third regime: 

e‘^ <C ||k'-||in <C 

in which 74 and A are negligible, and the Fermi surface is approximated by {p^j}je{o,i, 2 , 3 }) 
around which the Schwinger function diverges linearly. 

Remark: If 74 = AO, then the error term 0(e^||k) ||j^j) in (1.12) vanishes identically, 
which allows us to extend the third regime to all momenta satisfying 

||k'-||ni < 


1.3 Main Theorem 

We now state the Main Theorem, whose proof will occupy the rest of the paper. 
Roughly, our result is that as long as |t/| and e are small enough and 74 = A = 0 (see 
the remarks following the statement for an explanation of why this is assumed), the 
free energy and the two-point Schwinger function are well dehned in the thermodynamic 
and zero-temperature limit |A|,/3 —)• 00 , and that the two-point Schwinger function is 
analytically close to that with 1/ = 0. The effect of the interaction is shown to merely 
renormalize the constants of the non-interacting Schwinger function. 

We define 

B^:=Rx{R^/{ZG,+ZG 2 )), Gi:= ^2 := 

where the physical meaning of M^/(ZGi -|- ZG 2 ) is that of the first Brillouin zone, and 
G 17 are the generators of the dual lattice. 
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I-Main Theorem- 

If 74 = A = 0, then there exists Uq > 0 and eo > 0 such that for all \U\ < Uq and e < eo, 
the specific ground state energy 

eo ;= - lim lim ^ log(Tr(e“^^)) 

/3^oo |A|^oo p|A| 

exists and is analytic in U. In addition, there exist eight Fermi points j}aj=±j=o,i, 2,3 
such that: 

^^,0 = Pf,o> Ipfj - PfjI ^ (const.) \U\e^, j = 1, 2, 3, (1.13) 

and, Vk G Boo \ {p'pj}i.j=±,i=o,i, 2 , 3 , the thermodynamic and zero-temperature limit of 
the two-point Schwinger function, lim_s_^oo hm| 4 |_^oo 'S 2 (k), exists and is analytic in U. 


Remarks: 

• The theorem requires 74 = A = 0. As we saw above, those quantities play a 
negligible role in the non-interacting theory as long as we do not move beyond the 
third regime. This suggests that the theorem should hold with 74 , A 7 ^ 0 under 
the condition that /3 is not too large, i.e., smaller than (const.) e“^. However, that 
case presents a number of extra technical complications, which we will spare the 
reader. 

• The conditions that \U\ < Uq and e < eo are independent, in that we do not 
require any condition on the relative values of \U\ and e. Such a result calls for 
tight bounds on the integration over the first regime. If we were to assume that 
\U\ <C e, then the discussion would be greatly simplified, but such a condition 
would be artificial, and we will not require it be satisfied. L. Lu [Lul3] sketched 
the proof of a result similar to our Main Theorem, without discussing the first two 
regimes, which requires such an artificial condition on U/e. The renormalization 
of the secondary Fermi points is also ignored in that reference. 


In addition to the Main Theorem, we will prove that the dominating part of the 
two point Schwinger function is qualitatively the same as the non-interacting one, with 
renormalized constants. This result is detailed in Theorems l.I, 1.2 and 1.3 below, each 
of which refers to one of the three regimes. 

1 - First regime. Theorem 1.1 states that in the first regime, the two-point 
Schwinger function behaves at dominant order like the non-interacting one with renor¬ 
malized factors. 
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I-Theorem 1.1- 

Under the assumptions of the Main Theorem, if Ce ^ ||k — p'polli ^ ® suitable 

C > 0, then, in the thermodynamic and zero-temperature limit. 


S 2 (k) = 


koko + kP 


/ -iko o_ q e \ 

0 —iko ^ 0 

0 r -ih 0 

V ^ 0 0 -iko I 


(l+r(k)) (1.14) 


where 


+ k') = 0((l |U|| log llk'llil) llk'lli, e||k'||i), 

and, for (/cq, k'^, k'y) := k - 

- 3 

ko:=ziko, ko:=ziko, ( :=-viiik'^ + uky) 

in which {zi,zi,vi) G satisfy 

|l-zi| ^Ci|U|, |l-zi| ^Ci|U|, |l-ui| ^Ci|C/| 


(1.15) 


(1.16) 


(1.17) 


for some constant Ci > 0 (independent of U and e). 


Remarks: 

• The singularities of S 2 are approached linearly in this regime. 

• By comparing (1.14) with its non-interacting counterpart (3.8), we see that the 
effect of the interaction is to renormalize the constants in front of ko and ^ in (3.8). 

• The inter-layer correlations, that is the {a, b} x {a, b} components of the dominat¬ 
ing part of S 2 (k) vanish. In this regime, the Schwinger function of bilayer graphene 
behave like that of two independent graphene layers. 


2 - Second regime. Theorem 1.2 states a similar result for the second regime. 
As was mentioned earlier, two of the components are massive in the second (and third) 
regime, and we first perform a change of variables to isolate them, and state the result 
on the massive and massless components, which are denoted below by and 
respectively. 
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I-Theorem 1.2- 

Under the assumptions of the Main Theorem, if ^ ||k —p'^qUii ^ C~^e for a suitable 
C > 0, then, in the thermodynamic and zero-temperature limit, 


S2(k) 


1 M(k)t 
0 1 



0 


1 0 
M(k) 1 


(1 -4 r(k)) 


where: 




4”’(k) = 


'i-liko 7i(C 

f 2 


II 

* \2 


V 7ir iliko 




) •’2 


S'r, - - 


-1 


0 

7f^ 0 


1 


and, for {ko, k'^, k'y) := k - p'^Q, 


7i := hT.27i, ^0 := Z2kQ, ^ := -V2{ik'^ + ojk'y) 

in which {rh 2 , Z 2 ,V 2 ) G satisfy 

|l-m 2 | ^ C2IUI, II-Z2I ^ C2IUI, II-U2I ^ C2IC/I 
for some constant 6*2 > 0 (independent of U and e). 


(1.18) 

(1.19) 

( 1 . 20 ) 

( 1 . 21 ) 

( 1 . 22 ) 

(1.23) 


Remarks: 

• The massless components {d, b} are left invariant under the change of basis that 
block-diagonalizes S 2 . Furthermore, M is small in the second regime, which implies 
that the massive components are approximately {a,b}. 

• As can be seen from (1.20), the massive part s^'^ of S 2 is not singular in the 

neighborhood of the Fermi points, whereas the massless one, i.e. is. 

• The massless components of S 2 approach the singularity quadratically in the spatial 
components of k. 

• Similarly to the first regime, by comparing (1.20) with (3.18), we find that the 
effect of the interaction is to renormalize constant factors. 


3 - Third regime. Theorem 1.3 states a similar result as Theorem 1.2 for the 
third regime, though the discussion is made more involved by the presence of the extra 
Fermi points. 
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I-Theorem 1.3- 

For j = 0,1, under the assumptions of the Main Theorem, if ||k — p'^jHin ^ C~^e^ for 
a suitable C > 0 , then 


S 2 (k) = 


1 M(k)t 

0 1 


,(M) 


z;(™) 


(k) 


M(k) 1 


(1 + r(k)) 


where 


r{p%,j + k'-) = 0(e + e| log ||k'Jiii||17|), e(l + | loge||[/|)) 


_(m) 


(k) = 


iko,j 73X* 

^o,i+ 73l^iP V 73% ikoj 

M(k) := -3^ 
71,j 


- _ 
*0 — 


0 7i:i 

7i,l 0 


H* 0 

j _ 

0 


and, for (%, k'y) := k - 

%,i := Z3,jko, 71,i := 1 ^ 3 ^ 71 , xq := i;3,o|(i% - u)ky) =: -E^ 

xi := I ( 3 u 3 ,iiA:^ + ws^iujk'^) , Si := 111,3,17173 + n^iik'^ + Ws^iky 
in which Z 3 j,V 3 j,V 3 j,W 3 j,W 3 j) € satisfy 

l™'3,i - 1| + 1^-37 - 1| ^ C3\U\, \Z3J - 1| ^ <731171, 

|73,i - 1| + \v3,j - 1| ^ <73|17|, \w3j - 1| + |lh3,i - 1| ^ C3\U\ 
for some constant 6*3 > 0 (independent of U and e). 


(1.24) 

(1.25) 

(1.26) 

(1.27) 

(1.28) 

(1.29) 


Theorem 1.3 can be extended to the neighborhoods of p'^^ with j = 2,3, by taking 
advantage of the symmetry of the system under rotations of angle 2TTj3: 

I-Extension to 7 = 2,3- 

For j = 2,3, under the assumptions of the Main Theorem, if ||k — p'^jHin ^ C~^e^ for 
a suitable C > 0 , then 


S 2 (k'' + p^o) — 



S2(Tk) + p*^ 


J-U)j 


0 % 


Tk'+p-^_ 


(1.30) 


where T(/co, kx, ky) denotes the rotation of the kx and ky components by an angle 27r/3, 
T(ko,k^,ky) ■= and p ^_4 = p^^. 


Remarks: 

• The remarks below Theorem 1.2 regarding the massive and massless fields hold 
here as well. 
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• The massless components of S 2 approach the singularities linearly. 

• By comparing (1.24) with (3.25) and (3.32), we find that the effect of the interaction 
is to renormalize the constant factors. 


1.4 Sketch of the proof 

In this section, we give a short account of the main ideas behind the proof of the 
Main Theorem. 

1 - Multiscale decomposition. The proof relies on a multiseale analysis of 
the model, in which the free energy and Schwinger function are expressed as successive 
integrations over individual scales. Each scale is defined as a set of k’s contained inside 
an annulus at a distance of 2^ for h G Z around the singularities located at p'pj- The 
positive scales correspond to the ultraviolet regime, which we analyze in a multiscale 
fashion because of the (very mild) singularity of the free propagator at equal imaginary 
times. It may be possible to avoid the decomposition by employing ideas in the spirit 
of [PS08]. The negative scales are treated differently, depending on the regimes they 
belong to (see below), and they contain the essential difficulties of the problem, whose 
nature is intrinsically infrared. 

2 - First regime. In the first regime, i.e. for —I ^ h ^ he ;= log 2 e, the 
system behaves like two uncoupled graphene layers, so the analysis carried out in [GMIO] 
holds. From a renormalization group perspective, this regime is super-renormalizahle: 
the scaling dimension of diagrams with 21 external legs is 3 — 21 , so that only the two- 
legged diagrams are relevant whereas all of the others are irrelevant (see section 5.2 for 
precise definitions of scaling dimensions, relevance and irrelevance). This allows us to 
compute a strong bound on four-legged contributions: 

|ITf^(k)| ^ (const.) \U\2‘^’^ 

whereas a naive power counting argument would give \U\ 2 ^ (recall that with our con¬ 
ventions h is negative). 

The super-renormalizability in the first regime stems from the fact that the Fermi 
surface is 0-dimensional and that Hq is linear around the Fermi points. While performing 
the multiscale integration, we deal with the two-legged terms by incorporating them 
into Hq, and one must therefore prove that by doing so, the Fermi surface remains 0- 
dimensional and that the singularity remains linear. This is guaranteed by a symmetry 
argument, which in particular shows the invariance of the Fermi surface. 

3 - Second regime. In the second regime, i.e. for 3/ie <C h <C the singu¬ 
larities of Hq are quadratic around the Fermi points, which changes the power eounting 
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of the renormalization group analysis; the scaling dimension of 2Z-legged diagrams be¬ 
comes 2 — I so that the two-legged diagrams are still relevant, but the four-legged ones 
become marginal. One can then check [ValO] that they are actually marginally relevant, 
which means that their contribution increases proportionally to \h\. This turns out not 
to matter; since the second regime is only valid for h ^ may only increase by 

31/if I, and since the theory is super-renormalizable in the first regime, there is an extra 
factor 2^' in so that actually increases from 2^' to 3|/i£|2^% that is to say 

it barely increases at all if e is small enough. 

Once this essential fact has been taken into account, the renormalization group 
analysis can be carried out without major difficulties. As in the first regime, the invari¬ 
ance of the Fermi surface is guaranteed by a symmetry argument. 

4 - Third regime. In the third regime, i.e. for h <C 3/ie, the theory is again 
super-renormalizable (the scaling dimension is 3—21). There is however an extra difficulty 
with respect to the first regime, in that the Fermi surface is no longer invariant under 
the renormalization group flow, but one can show that it does remain 0-dimensional, and 
that the only effect of the multiscale integration is to move along the line between 
itself and PpQ- 


1.5 Outline 

The rest of this paper is devoted to the proof of the Main Theorem and of Theo¬ 
rems 1.1, 1.2 and 1.3. The sections are organized as follows. 

• In section 2, we define the model in a more explicit way than what has been done 
so far; then we show how to compute the free energy and Schwinger function 
using a Fermionic path integral formulation and a determinant expansion, due to 
Battle, Brydges and Federbush [BF78, BF84], see also [BK87, AR98]; and finally 
we discuss the symmetries of the system. 

• In section 3, we discuss the non-interacting system. In particular, we derive de¬ 
tailed formulae for the Fermi points and for the asymptotic behavior of the prop¬ 
agator around its singularities. 

• In section 4, we describe the multiscale decomposition used to compute the free 
energy and Schwinger function. 

• In section 5, we state and prove a power counting lemma, which will allow us to 
compute bounds for the effective potential in each regime. The lemma is based on 
the Gallavotti-Nicolo tree expansion [GN85], and follows [BG90, GMOl, GilO]. We 
conclude this section by showing how to compute the two-point Schwinger function 
from the effective potentials. 

• In section 6, we discuss the integration over the ultraviolet regime, i.e. scales h > 0. 
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• In sections 7, 8 and 9, we discuss the multiscale integration in the first, second 
and third regimes, and complete the proofs of the Main Theorem, as well as of 
Theorems 1.1, 1.2, 1.3. 


2 The model 


From this point on, we set 74 = AO. 

In this section, we define the model in precise terms, re-express the free energy 
and two-point Schwinger function in terms of Grassmann integrals and truncated expec¬ 
tations, which we will subsequently explain how to compute, and discuss the symmetries 
of the model and their representation in this formalism. 


2.1 Precise definition of the model 

In the following, some of the formulae are repetitions of earlier ones, which are 
recalled for ease of reference. This section complements section 1.1, where the same defi¬ 
nitions were anticipated in a less verbose form. The main novelty lies in the momentum- 
real space correspondence, which is made explicit. 


1 - Lattice. As mentioned in section 1, the atomic structure of bilayer graphene 
consists in two honeycomb lattices in so-called Bernal or AB stacking, as was shown 
in figure 1.1. It can be constructed by copying an elementary cell at every integer 
combination of 


3 v/3 


v/3 




( 2 . 1 ) 


where we have chosen the unit length to be equal to the distance between two nearest 
neighbors in a layer (see figure 2.1). The elementary cell consists of four atoms at the 
following coordinates 

(0,0,0); (0,0,c); (-l,0,c); (1,0,0) 


given relatively to the center of the cell, c is the spacing between layers; it can be 
measured experimentally, and has a value of approximately 2.4 [TMe92]. 


We define the lattice 

A := [mh + n 2 l 2 , (ni, 712 ) G {0, • • • , L - 1}^} (2.2) 

where L is a positive integer that determines the size of the crystal, that we will even¬ 
tually send to infinity, with periodic boundary conditions. We introduce the intra-layer 
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Figure 2.1: decomposition of the crystal into elementary cells, represented by the blue 
(color online) rhombi. There are four atoms in each elementary cell: 9 of type a at 
(0,0,0), O of type b at (0, 0, c), Q of type a at (—1,0, c) and ■ of type b at (0,0, c). 


nearest neighbor vectors: 


:= ( 1 , 0 , 0 ), 



53 := 




(2.3) 


The dual of A is 

(2.4) 

(2.5) 

It is defined in such a way that Vx G A, V/c G A, 

^ikxL _ 


A := + ^^ 2 , {mi,m 2 ) G {0, 


,L- 


ir} 


with periodic boundary conditions, where 


27r 27r 


3 ’ V3’ 


/ 27r 27r 


Cl — —, ^,0 , G2 — { -7^, 0 


V 3 ’ V3’ 


Since the third component of vectors in A is always 0, we shall drop it and write vectors 
of A as elements of M^. In the limit L —)• 00 , the set A tends to the torus Aqo = 
M?/{ZGi +ZG 2 ), also called the Brillouin zone. 
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2 - Hamiltonian. Given x G A, we denote the Fermionic annihilation operators 
at atoms of type a, b, a and b within the elementary cell centered at x respectively by ax, 
bx, dx-Si and bx+Si- The corresponding creation operators are their adjoint operators. 

We recall the Hamiltonian (1.1) 

n = no + ni 

where T-Lq is the free Hamiltonian and Hi is the interaction Hamiltonian. 


2-1 - Free Hamiltonian. As was mentioned in section 1, the free Hamilto¬ 
nian describes the hopping of electrons from one atom to another. Here, we only consider 
the hoppings 7o, 7i, 73 , see figure 1 . 2 , so that Hq has the following expression in x space: 

Hq := -70 (albx+Sj + -t bldxSj + - 7i X] 

oiGA xGA 

J=l,2,3 

xGA 

i=l,2,3 

( 2 . 6 ) 

Equation (2.6) can be rewritten in Fourier space as follows. We define the Fourier 
transform of the annihilation operators as 

hfc := Y1 ’ ~^k-=Yl , ak:=Y^ e^^^dxSi > bk:=Y^ e^^^bx+s^ 

xGA xGA fcGA xGA 

(2.7) 

in terms of which 

= ( 2 . 8 ) 

fceA 

where |A| = L^, Ak is a column vector, whose transpose is = {dk,bk,dk,bk), 


Hoik) := 


0 7i 0 7on*(A:) 

7i 0 7o^ik) 0 

0 -fon*ik) 0 73ll(A:)e3**^- 

7ofl(A:) 0 73ll*(A:)e-3*^- 0 


(2.9) 


n{k) := = 1 + 2 e-*i^ 


We pick the energy unit in such a way that 70 = 1. 



2-2 - Interaction. We now define the interaction Hamiltonian. We first 
define the number operators nf for a G {a, b, d, b} and x G A in the following way: 


< = alox , 


G ~ blfix , n^ — d^^_^^dx-Si , — bl,_^g_^bx+5i 


( 2 . 10 ) 
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and postulate the form of the interaction to be of an extended Hubbard form: 

Hi — U v{x + da - y - da') (n^'- (2.11) 

(a;,j/)eA2 {a,a')£{a,b,a,bP ^ ^ ^ ^ 

where the da are the vectors that give the position of each atom type with respect to 
the centers of the lattice A: da ■= 0, d^ := (0,0, c), da ■= (0,0, c) — (5i, db := (5i and v 
is a bounded, rotationally invariant function, which decays exponentially fast to zero at 
inhnity. In our spin-less case, we can assume without loss of generality that u(0) = 0. 


2.2 Schwinger function as Grassmann integrals and expectations 


The aim of the present work is to compute the specific free energy and the two-point 
Schwinger function. These quantities are dehned for finite lattices by 

1 


/a := -T^log Tr e 


m 

where /3 is inverse temperature and 

'Sa',a(xi - X2) := (T(a( 


-m 


Xi ^X2 ^ 


Tr(e-/^^T«4,)) 

Tr(e-/3«) 


( 2 . 12 ) 


(2.13) 


in which (a, a') G := {a, 6, a, 6}^; xi ^2 = (h, 2 , 3 ^ 1 , 2 ) with ti ^2 G [0,/3); ax = 
and T is the Fermionic time ordering operator defined in (1.8). Our strategy 
essentially consists in deriving convergent expansions for f\ and s, uniformly in |A| and 
/?, and then to take /3, |A| —)■ 00 . 


However, the quantities on the right side of (2.12) and (2.13) are somewhat diffi¬ 
cult to manipulate. In this section, we will re-express /a and s in terms of Grassmann 
integrals and expectations, and show how such quantities can be computed using a de¬ 
terminant expansion. This formalism will lay the groundwork for the procedure which 
will be used in the following to express f\ and s as series, and subsequently prove their 
convergence. 


1 - Grassmann integral formulation. We first describe how to express (2.12) 
and (2.13) as Grassmann integrals. The procedure is well known and details can be found 
in many references, see e.g. [GMIO, appendix B] and [GilO] for a discussion adapted 
to the case of graphene, or [GMOl] for a discussion adapted to general low-dimensional 
Fermi systems, or [BG95] and [Sall3] and references therein for an even more general 
picture. 


1-1 - Definition. We first define a Grassmann algebra and an integration 
procedure on it. We move to Fourier space: for every a G A := {a, b, d, b}, the operator 
a(^t,x) is associated 


Ak=(A:o,fc) 


P Jo 
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with /cq G 27r/3“^(Z + 1/2) (notice that because of the 1/2 term, feo / 0 for finite /3). We 
notice that k G fi/ 3 ,L := (27r/3“^(Z + 1/2)) x A varies in an infinite set. Since this will 
cause trouble when defining Grassmann integrals, we shall impose a cutoff M G N: let 
Xo(/o) be a smooth compact support function that returns 1 if p ^ 1/3 and 0 if p ^ 2/3, 
and let 

BIl ■= n {(fco, k), xo(2-'^|A:o| / 0)}. 

To every (dkjdj^) for a G ^ and k G .B/i we associate a pair of Grassmann variables 
(^ko’^ica)’ consider the finite Grassmann algebra (i.e. an algebra in which the 

■0 anti-commute with each other) generated by the collection {V’kalkle* ■ define 
the Grassmann integral 

/ aeA 

as the linear operator on the Grassmann algebra whose action on a monomial in the 
variables 0^^, is 0 except if said monomial is OkeB^ ^ ^ka'*^ka ^ permutation of 

the variables, in which case the value of the integral is determined using 


aeA 


a&A 


n n '^k,a^k,, 


keZ3 


/3.L 


,keB 


/3,L 


= 1 


(2.14) 


along with the anti-commutation of the 0. 


In the following, we will express the free energy and Schwinger function as Grass¬ 
mann integrals, specified by a propagator and a potential. The propagator is a 4 x 4 
complex matrix p(k), supported on some set B C B^ and is associated with the Gaus¬ 
sian Grassmann integration measure 

Pgi.di’) ■■= ^n(/3|A|)4detp(k) c?0k,a#k,a^^ exp ^kFHk)00 

(2.15) 

Gaussian Grassmann integrals satisfy the following addition prineiple: given two prop¬ 
agators gi and g 2 , and any polynomial ip(0) in the Grassmann variables, 

I Pg,+g,{d^) ^(0) = I PgMi^l) j ^§,(^02) ^1(01 +02)- (2.16) 



1-2 - Free energy. We now express the free energy as a Grassmann integral. 
We define the free propagator 

p^M(k) := xoi2-^\ko\){-ikol - Hoik))-^ (2.17) 
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and the Gaussian integration measure P<^M{dip) = Pg^j^{dip). One can prove (see e.g. 
[GMIO, appendix B]) that if 

^log J ( 2 . 18 ) 

is analytic in U, uniformly as M —)■ oo, a fact we will check a posteriori, then the finite 
volume free energy can be written as 

/a = /o,A - Jim^ ^ log J P^Midi^) ( 2 . 19 ) 

where /o,a is the free energy in the ?7 = 0 case and, using the symbol f dx as a shorthand 
for 


V{i;) = U Y1 

{a,a')GA'^ 


dxdy Wa,a'(x - 


( 2 . 20 ) 


in which Wa,a'{^) '■= 
delta function, and 


(5(xo)u(x + da — da'), where S{xo) denotes the /^-periodic Dirac 


V’x.a := 




m 


<c 


^diikx 




( 2 . 21 ) 


Notice that the expression of V(V^) in (2.20) is very similar to that of Pj, with an added 
imaginary time (xo,yo) and the Ox replaced by 'ipx,a, except that (aiox — 1/2) becomes 
V’x aV'x o- Roughly, the reason why we “drop the 1/2” is because of the difference between 

the anti-commutation rules of and V’x,a (i-e., {ax,ax} = 1, vs. {V’xa^V'xol = 0)- 

More precisely, taking x = (xq, x) with xq G (—/3, /3), it is easy to check that the limit as 
M —> oo of (/,^m(x) := J P^M{d'ip)'^Px'^o is equal to s(x), if x / 0, and equal to s(0)-|-l/2, 
otherwise. This extra +1/2 accounts for the “dropping of the 1/2” mentioned above. 


1-3 - Two-point Schwinger fnnction. The two-point Schwinger function 
can be expressed as a Grassmann integral as well; under the condition that 


fP^uidi^) 


Oil 


f RgM(#) 


( 2 . 22 ) 


is analytic in U uniformly in M, a fact we will also check a posteriori, then one can prove 
(see e.g. [GMIO, appendix B]) that the two-point Schwinger function can be written as 


Sai,Q: 2 (k) = lim 

M—^OD 


fP^M(d^) e 


(2.23) 


/ P^M(dtp) 

In order to facilitate the computation of the right side of (2.23), we will first rewrite it 


as 


^ai^a2 (k) = lim / dX a dJ+ log / P^Mid^e 


+ 4,C2 


M^oo 


( 2 . 24 ) 
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where ^ and are extra Grassmann variables introduced for the purpose of the 

computation (note here that the Grassmann integral over the variables 

as a functional derivative with respect to the same variables, due to the Grassmann 

integration/derivation rules). We define the generating functional 

Wii.,AJ := V(^) - (2.25) 


2 - Expectations. We have seen that the free energy and Schwinger function 
can be computed as Grassmann integrals, it remains to see how one computes such 
integrals. We can write (2.18) as 


log / - ,V; =: E S^^lM(V;iV). 


(2.26) 


where the truncated expectation is dehned as 

gN 


gn p 

,V7v) := J P^M{di^) 


. (2.27) 


Ai=***=Ajv=0 


in which (Vi, • • • , Vat) is a collection of commuting polynomials and the index refers 
to the propagator of P^M{d'4’)- A similar formula holds for (2.22). 


The purpose of this rewriting is that we can compute truncated expectations in 
terms of a determinant expansion, also known as the Battle-Brydges-Federbush for¬ 
mula [BF78, BF84], which expresses it as the determinant of a Gram matrix. The 
advantage of this writing is that, provided we first re-express the propagator 5,gM(k) in 
x-space, the afore-mentioned Gram matrix can be bounded effectively (see section 5.2). 
We therefore first define an x-space representation for 5 (k): 


5=sm(x) 


1 

M 


E 


e**^'"ff^M(k). 


(2.28) 


The determinant expansion is given in the following lemma, the proof of which can be 
found in [GMOl, appendix A.3.2], [GilO, appendix B]. 


I-Lemma 2.1- 

Gonsider a family of sets P = (Pi, • • • , P^) where every Pj is an ordered collection 
of Grassmann variables, we denote the product of the elements in Pj by p. := Oi/jePj 

We call a pair (V’x V’x/ ^ ^ line, and define the set of spanning trees T(P) 

as the set of collections T of lines that are such that upon drawing a vertex for each 
Pi in P and a line between the vertices corresponding to Pi and to Pj for each line 
(V'x CO V'x' q') ^ P is such that '0x« ^ Pi V’x'a' ^ Pp ^kie resulting graph is a 
tree that connects all of the vertices. 
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For every spanning tree T G T(P), to each line I = (V’xa; V’x'«') ^ ^ we assign a 
propagator gi := ga,a'{^ - x'). 

If P contains 2(n + s — 1) Grassmann variables, with n G N, then there exists a 
probability measure dPT{t) on the set of n x n matrices of the form t = M'^M with M 
being a matrix whose columns are unit vectors of W^, such that 

dPrit) detG(^)(t) (2.29) 

reT(P) zeT 

where ctt G { — 1,1} and is an n x n complex matrix each of whose components 

is indexed by a line I 0 T and is given by 

= tigi 

(if s = 1, then T(P) is empty and both the sum over T and the factor ot OzeT 9i should 
be dropped from the right side of (2.29)). 


Lemma 2.1 gives us a formal way of computing the right side of (2.26). However, 
proving that this formal expression is correct, in the sense that it is not divergent, will 
require a control over the quantities involved in the right side of (2.29), namely the prop¬ 
agator g^M- Since, as was discussed in the introduction, g^M is singular, controlling the 
right side of (2.26) is a non-trivial task that will require a multiscale analysis described 
in section 4. 


2.3 Symmetries of the system 

In the following, we will rely heavily on the symmetries of the system, whose 
representation in terms of Grassmann variables is now discussed. 

A symmetry of the system is a map that leaves both 

ho :='^'iptg~^{^-y)'ipy (2.30) 

x,y 


and V{ip) invariant was defined in (2.20)). We define 



’(2.31) 


as well as the Pauli matrices 
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We now enumerate the symmetries of the system, and postpone their proofs to ap¬ 
pendix E. 


1 - Global U{1). 


For 6 G M/(27rZ), the map 




, ^±iel± 


(2.32) 


is a symmetry. 

2 - 27r/3 rotation. Let Tk := {ko,e~'^^^‘^k), I 2 := (3/2,—\/3/2) and 7k := 
g-i(i 2 -fc)o- 3 ^ the mapping 


^k 


^Tk 

• 'Trk'^Tk’ 


X+ 7"t 

‘^Tk 'Tk 


(2.33) 


is a symmetry. 


3 - Complex conjugation. The map in which 


and every complex coefficient of 
symmetry. 


<^k 


■<^-k- 


(2.34) 


ho and V is mapped to its complex conjugate is a 


4 - Vertical reflection. Let i?^k = {ko,ki, —k 2 ), 



is a symmetry. 


(2.35) 


5 - Horizontal reflection. Let 7?/ik = (/cq, —A:i, ^ 2 ), 

f ^k ' ^ ^k ' ^ ^Thk*^! 

1 '^k ' ^ ‘^k ' ^ ‘^Thk'^1 

is a symmetry. 

6 - Parity. Let Pk = {ko, —ki, — ^ 2 ), 


(2.36) 



is a symmetry. 


7 - Time inversion. Let Ik = (— /cQ) ^2); the mapping 

1 '^k ' ^ “^3<^7k’ '^ic ' ^ '^/k^3 

is a symmetry. 


(2.38) 


3 Free propagator 

In section 2.2, we showed how to express the free energy and the two-point 
Schwinger function as a formal series of truncated expectations (2.26). Controlling 
the convergence of this series is made difficult by the fact that the propagator g^M is 
singular, and will require a finer analysis. In this section, we discuss which are the sin¬ 
gularities of g^M and how it behaves close to them, and identify three regimes in which 
the propagator behaves differently. 


3.1 Fermi points 


The free propagator is singular if Iq = 0 and k is such that HQ{k) is not invertible. 
The set of such /c’s is called the Fermi surface. In this subsection, we study the properties 
of this set. We recall the definition of Hq in (2.9), 


Hoik) := 


/ 0 71 

7 i 0 

0 Q*ik) 

V nik) 0 


0 

Q{k) 

0 


n*ik) \ 

0 

73Ll(A:)e3*^- 

0 


so that, using corollary B.2 (see appendix B), 


det Hoik) 


n\k)-ji^3n*ik)e-^^^^ 


(3.1) 


It is then straightforward to compute the solutions of det Hoik) = 0 (see appendix A for 
details): we find that as long as 0 < 7173 < 2, there are 8 Fermi points: 


^, (u arccos 


1-7173 

2 


Pf,i 

Pf,2 - (f + i arccos (^ ^1+7173^-7173 ) 


„a;_ (277 


Pp^ ■= I — 3 arccos 


/ 71+7173(2-7173) 

I 2 


, arccos 


, arccos 


(l±m,)) 


( 1+7173 

I 2 


(3.2) 
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for u) G {—, +}. Note that 


Pf,i = Pf,o + + C>(e^), + (^^TiTs, -^^7173) + O(e^), 

Pf ,3 = PFfl + (-757173, -W57173) + 0 {e^). 

The points p‘^j for j = 1, 2, 3 are labeled as per figure 1.3. 


(3.3) 


3.2 Behavior around the Fermi points 

In this section, we compute the dominating behavior of 5 (k) close to its singu¬ 
larities, that is close to ppj := {0,p%j)- We recall that ^(k) ;= {—ikot + Ho{k)) and 
m = Xo{2-^\ko\)A-Hk). 

1 - First regime. We define k' := k — = (/c(,, ky), k' := {ko, k'). We have 

^(Pf,o + k') = -{ik!^ + ojk'y) + 0{\k'\^) =; C + 0{\k'\^) (3.4) 

so that, by using (B.2) with (a, b, c,}:,^) = —(71, f2(A:),73n(/c)e^*^"=, /cq, ko), 

det ^(p'^^Q -|- k') = (/cq -|- I^P)^ -|- 0(||k'||j, e^||k'||j) (3-5) 

where 

llk'lli ;= \lkl + |^|2 (3.6) 

in which the label i stands for “first regime”. If 

Kie < ||k'||i < kq (3.7) 

for suitable constants ki,ko > 0, then the remainder term in (3.5) is smaller than the 
explicit term, so that (3.5) is adequate in this regime, which we call the “first regime”. 
We now compute the dominating part of A~^ in this regime. The computation is carried 
out in the following way: we neglect terms of order 71 , 73 and \k'\‘^ in A, invert the 
resulting matrix using (B.3), prove that this inverse is bounded by (const.) ||k'||j^^, and 


deduce a bound 

on the error terms. We thus find 





0 

* 

1 

0 

0 

r \ 


1 

0 

-b 

1 

0 

-iko 


0 

(l + 0(||k'||i,e||k'||fi)) 

fc§ + |{P 

0 

r 

-iko 

0 


V C 

0 

0 

-iko ) 



and 

l^■^(PF,o + k')l ^ (const.) Ilk'llf^ (3.9) 
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Note that, recalling that the basis in which we wrote ^ is {a,b,d,b}, each graphene 
layer is decoupled from the other in the dominating part of (3.8). 

2 - Ultraviolet regime. The regime in which ||k'||i ^ kq for both a; = ±, and 
is called the ultraviolet regime. For such k' =: k — easily checks that 

|A“^(k)| ^ (const.) |k|~^. (3.10) 


3 - Second regime. We now go beyond the first regime: we assume that 
k'lli ^ /tie and, using again (3.4) and (B.2), we write 

det ^(p'^^Q + k') = 7 i/co + + 0(e'^/^||k^||j/ , e^||k'||n, e||k'||ji) (3.11) 


where 



If 

K 2 e^ < ||k'||n < ki e 


(3.12) 

(3.13) 


for suitable constants /t 2 , /ti > 0, then the remainder in (3.11) is smaller than the explicit 
term, and we thus define the “second regime”, for which (3.11) is appropriate. 


We now compute the dominating part of A~^ in this regime. To that end, we 
define the dominating part £,iiA of A by neglecting the terms of order 73 and jfc'p in A 
as well as the elements Aaa and Aj^^ (which are both equal to —iko), block-diagonalize 
it using proposition C.l (see appendix C) and invert it: 


£iii(k)) 


-1 


1 Mn(k)t 
0 1 


AM) 


^ii 


0 

0 a<7(k) 


1 0 
Mn(k) 1 


where 


A^l — 


qi 


0 7]) 

7f^ 0 


-1 


®ii ^(Pf,o + 1^0 •— 


71 


and 


M„(p"p,. + k') :=-!(«• " 


We then bound the right side of (3.14), and find 
I (^£iii(p)(o + k')) I ^ (const.) 


niko (? 


* '\2 


7?fco + 1^1"^ V U hiko 


-1 
e-V2||k 


e " e 

/11-1/2 

II 


-V2||U||-V2 

iik'iir/ 


(3.14) 


(3.15) 


(3.16) 


(3.17) 
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in which the bound should be understood as follows; the upper-left element in (3.17) is 
the bound on the upper-left 2x2 block of (i2iiA)“^, and similarly for the upper-right, 
lower-left and lower-right. Using this bound in 


i-^(k) = (i2iii(k)) ' [l + (i(k) - £iii(k)) (£iii(k)) 


-ii-i 


we deduce a bound on the error term in square brackets and find 


"4 (p‘f,o + k') — 


1 Mn(k)t ( ajf) 0 


0 


(m) 


which implies the analogue of (3.17) for A 
I^“^(Pf,o+ k')l ^ (const.) { 


)/ 

' 1 

0 ^ 


(k) j 1 

, Mii(k) 

1 ) 


•(1 + 0(e-i/2||k/ 






(3.18) 

-1 

711-1/2 

^ llll 

e-i/2||k'|| 

Ilk'IlH 

-1/2 \ 

" )■ 

(3.19) 


Remark: Using the explicit expression for A ^(p'pQ + k') obtained by applying propo¬ 
sition B.l (see appendix B), one can show that the error term on the right side of (3.18) 
can be improved to 0(e“^||k'||ii, e^/^||k'||jj^'^^)). Since we will not need this improved 
bound in the following, we do not belabor further details. 

4 - Intermediate regime. In order to derive (3.18), we assumed that llk'llii < 
Rie with Ri small enough. In the intermediate regime defined by Rie < ||k'||ii and 
||k'||i < Kie, we have that ||k'||i ~ ||k'||ii ~ e (given two positive functions a(e) and 6(e), 
the symbol a ~ 6 stands for cb ^ a ^ Cb for some universal constants C > c > 0). 
Moreover, 

det yl(p‘^g -|- k') = (fcg -|- |^P)^ -|- 7i fcg + 0{e^) (3.20) 

therefore | det^l > (const.) e^ and 

l^~^(p‘F,o>k')| ^ (const.) e~^ (3.21) 

which is identical to the bound at the end of the first regime and at the beginning of the 
second. 

5 - Third regime. We now probe deeper, beyond the second regime, and assume 
that ||k'||n ^ K 2 (^- Since we will now investigate the regime in which |/c'| < (const.) e^, 
we will need to consider all the Fermi points p^j with j g{0,1,2,3}. 

5-1 - Around p'pg. We start with the neighborhood of p'pg: 

det A(p‘p^o + kO = 7i (^0 + tII^P) + 0(e ^||k'||in) (3.22) 
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where 


(3.23) 


llk'lliii := \Jkl + 73ICP- 

The third regime around p'^Q is defined by 

||k^||iii < (3.24) 

for some K 2 < K 2 - The computation of the dominating part of A~^ in this regime around 
Pfo is similar to that in the second regime, but for the fact that we only neglect the 
terms of order \k'\‘^ in A as well as the elements Aaa and Aj^^. In addition, the terms 
that are of order e“^||k'||jjj that come out of the computation of the dominating part of 
A in block-diagonal form are also put into the error term. We thus find 


i-i(k) = 


where 

- 

“■111,0 

and 


and 


1 Tfni,o(k)''’ 
0 1 


AM) 

hii,o 


0 i V Tfiii,o(k) 1 


1 0 


(l+0(e-=||k'||m)) 


-1 


0 

0 


)’ “m,o(PF,o + k'):- ^2+^2|^|2( 


Tfiii,o(PF,o + ( 0 ^ 


-iko 73^ 
73^ -iko 


1^ nP‘F,o + k')l ^ (const.) ^ ^_2 ^ . 


(3.25) 

(3.26) 

(3.27) 

(3.28) 


5-2 - Around p'^^. We now discuss the neighborhood of v^i- We define 
k'l := k — p^-^ = {k'l^^, k'l y) and k'^ ;= (/cq, k'^). We have 

+ k'l) = 7173 -|- Cl + 0{^\k'i\) (3.29) 

where 

6 := \{ik'i^^ + ujk'i y). 

Using (B.2) and (B.4), we obtain 

deti(p)(,i + k'l) = ^Ikl + |02 _ 717357*6-3*^1.-12 + 0(e^|A:oP) (3.30) 

where is evaluated at -|- k'l- Injecting (3.29) into this equation, we find 

det A(p‘^^i -I- k'l) = 7 i(A:o + TsI^iP) + 0(e^||k'i||fjj, e ^||k'i||jji) (3.31) 
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where 


3 

Xi := -{m'^ ,. + ujk[ y). 

The third regime around ^ is therefore defined by 

111 / 11 - — 3 

||ki||iii < K2 e 


where K 2 can be assumed to be the same as in (3.24) without loss of generality. The 
dominating part of A~^ in this regime around is similar to that around except 
that we neglect the terms of order e^k'i in A as well as the elements Aaa and Ayj^. As 
around p'pg, the terms of order e“^||k'^||jjj are put into the error term. We thus find 


i-i(k) 


1 Afni,i(k)'l’ 
0 1 


where 


0 


0 



1 0 \ 

Af[ii,i(k) 1 ) 

•(1 + 0(e,€“^||k'||ni)) 
(3.32) 


- 


0 7r' 

7i ^ 0 


(m) 

%I,1 


(Pf,i + kO 


1 

^0 +7|kip 


iko 733^1 
73x1 iko 


(3.33) 


and 


and 


-^iii.i (Pf,i + k]^) 


-73I - 


1 

71 


et 0 \ 
0 a J 


1 ^ ^(Pf.i +ki)| ^ (const.) 


p2||k' 11-1 

<=llk' 11“^ 


Ilk'illEl J ■ 


(3.34) 

(3.35) 


5-3 - Around p‘pj. The behavior of g(k) around ppj for j € {2,3} can be 
deduced from (3.32) by using the symmetry (2.33) under 27r/3 rotations: if we define 
kj := k — = (fe' a;, k'j y), k'- ;= (feg, k'^) then, for j = 2,3 and u±, 




1 0 
0 


A (Tk) + PF,7-a 


1 0 

n 

^ 'Tk'+p“ 


(3.36) 


where T and 7k were defined above (2.33), and p^4 = Pp^^. In addition, if k '2 and kg 
are in the third regime, then + O(e^). 


6 - Intermediate regime. We are left with an intermediate regime between 
the second and third regimes, defined by 

< ||k'||iii , ||k'||n < ^26^ and it2e^ < ||k' ||ni, Vj G {1,2,3}, (3.37) 
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which implies 


Ik'l 


III 


Ik'l 


II 


Ik'l 


III ~ e 


and 

det + k') = -ilkl + 1^2 - 7173^1^ + 0(e^°)- (3.38) 

One can prove (see appendix D) that injecting (3.37) into (3.38) implies that | det A| ^ 
(const.) e®, which in turn implies that 


1 ^ ^(Pf,o + k')! ^ (const.) 




(3.39) 


which is identical to the bound at the end of the second regime and at the beginning of 
the third. 


7 - Summary. Let us briefly summarize this sub-section: we defined the norms 

III := \/^o + 7 |I^P, ( 3 . 40 ) 

and identified an ultraviolet regime and three infrared regimes in which the free propa¬ 
gator 5(k) behaves differently; 

• for ||k'||i > kq, ( 3 . 10 ) holds. 

• for Kie < ||k'||i < Ko, ( 3 . 8 ) holds. 

• for K2e^ < ||k'||n < Rie, ( 3 . 18 ) holds. 

• for ||k'||ni < ^2^^, ( 3 . 25 ) holds, for ||k)|| III < K2e^, ( 3 . 32 ) holds, and similarly for 
the j = 2,3 cases. 



4 Multiscale integration scheme 

In this section, we describe the scheme that will be followed in order to compute 
the right side of (2.26). We will first define a multiscale decomposition in each regime 
which will play an essential role in showing that the formal series in (2.26) converges. In 
doing so, we will define effective interactions and propagators, which will be defined in k- 
space, but since we wish to use the determinant expansion in lemma 2.1 to compute and 
bound the effective truncated expectations, we will have to define the effective quantities 
in x-space as well. Once this is done, we will write bounds for the propagator in terms 
of scales. 
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4.1 Multiscale decomposition 

We will now discuss the scheme we will follow to compute the Gaussian Grass- 
mann integrals in terms of which the free energy and two-point Schwinger function were 
expressed in (2.19) and (2.24). The main idea is to decompose them into scales, and 
compute them one scale at a time. The result of the integration over one scale will then 
be considered as an effective theory for the remaining ones. 

Throughout this section, we will use a smooth cutoff function Xo(p)) which returns 
1 for p ^ 1/3 and 0 for p ^ 2/3. 


1 - Ultraviolet regime. Let f)o := [log2(Ko)J (in which Rq is the constant that 
appeared after (3.40) which defines the inferior bound of the ultraviolet regime). For 
h G {1)0, • • • , M} and h' G {^o + 1, • • • , M}, we define 

/=£h'(k) := Xo(2“^'|/co|), /,g(jo(k) := Xo(2“'’“||k - p‘/._ol|i), 

A'(k) := /^/i'(k) - /^/i'-i(k) (4.1) 

■= n suppf^h, ■= n suppA/, 

in which || • ||i is the norm defined in (3.40). In addition, we define 

fffc'(k) := A'(k)i"^(k), g^h{^) ■= /s£ft(k)i"^(k) (4.2) 

so that, in particular, 

p,gM(k) = p,gM-i(k) -h PM(k). 

Furthermore, it follows from the addition property (2.16) that for all h G {f}o, • • • , M—1}, 
< (4.3) 

where P^hidip^^^^) = 

-/3|A|Ffc-VW(V’(<^)) := -/3\A\Fh+i+log j 


OO (-.-.N 

= -mFh+i + E + V’^^"^);A^) 

7 V =1 

(4.4) 
.nd 

-m{Fh - Fh+l) - 

■= E ^-4p4.+i)(w<'‘+‘>«><'‘+'> + ik,a);JV) 

N=1 

(4.5) 
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in which the induction is initialized by 


;= V, ;= W, Fm := 0. 


2 - First regime. We now decompose the first regime into scales. The main 
difference with the ultraviolet regime is that we incorporate the quadratic part of the 
effective potential into the propagator at each step of the multiscale integration. This is 
necessary to get satisfactory bounds later on. The propagator will therefore be changed, 
or dressed, inductively at every scale, as discussed below. 

Let 1)1 := |'log2(/iie)] (in which ki is the constant that appears after (3.40) which 
defines the inferior bound of the first regime), and || • ||i be the norm defined in (3.40). 
We define for h G {f)i, • • • , f)o} and h' G {f)i + 1, • • • , f}o}, 

Uh,W^) ■= Xo(2"'*||k - ■= Uh'A'^) - 

^ suppf,^h,io, ■= ^li,L n s.\xppf^h>,ui 

and 

gh',Lj{^) ■= 9^h,U^) ■= 

For h G {1)1, • • • , ^0 ~ 1}) we define 

-f3\A\{Fh - Fh+i) - 

1 1\7V 

E ^4p4Li(v<'*+‘'G<'*+‘' + >/-«'■');'V) 

iV=l 

g(fjo)(^(s£fjo)) _|_ 'p(fjo)('^(s£fio)^ y(iio)(-.0(=£fjo)^ 

(4.8) 

and 

-miFh - Fh+i) - ik,a) 

N=l 

Q(^o)(^(^ilo)) +yp(i)o)(^(<f)o)^ := Jk a) 

(4.9) 

in which is quadratic in the ■0; is at least quartic and has no terms that 
are both quadratic in 0 and constant in Jk,Q; and is the truncated expectation 

defined from the Gaussian measure _ ('^0^'*"^^); iii which gh+i,uj is 

the dressed propagator and is defined as follows. Let (k) denote the kernel of 

c''“'(V’<«>) =: E E i’LT'LiUW'C’L W'lo) 


(4.6) 

(4.7) 
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(remark: the index in „ is redundant since given k, it is defined as the unique 
io that is such that k G it will however be needed when defining the x-space 

counterpart of below). We define gh,uj and by induction: := 

and, for h G {f)i + 1,..., f)o}, 

\ ■■= [hhMT^ + "'^(k) 

in which /</J^(k) is equal to 1/f^h,ujiy) if /^/i,a;(k) / 0 and to 0 if not. 

The dressed propagator is thus defined so that 


M 


[di,) I 




(4.11) 


(4.12) 


in which P^h = P^^^ Equation (4.11) can be expanded into a 

more explicit form: for h' G {f)i + 1,..., f)o} and h G {f)i, • • • , ^o}; 

l/i',a;(k) = fh'A^) (^h',c^(k)) , Yh,uj{^) = Uh,0J (ift,a;(k)^ (4.13) 

where 

to 

Ah,M:=A(k) + Uh,MwY(M)+ Y1 (4-14) 

h'=h+l 

(in which the sum should be interpreted as zero if h = f)o)- 


3 - Intermediate regime. We briefly discuss the intermediate region between 
regimes 1 and 2. We define 

PuM :=yo(2“'’M|k-PF,olli) -Xo(2“^M|k-PF,ollii) =: Uhi,ujO^) - Ui,uY'^) (4-15) 

where f)i := [log2(«:ie)J, from which we define 5[,^^^(k) and 5,g(j^^^(k) in the same way 
as in (4.13) with 

to 

ijj^^Jk) := i(k) + /^i^^,Jk)W»^)(k) + TEf ^(k). (4.16) 

The analogue of (4.12) holds here as well. 
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4 - Second regime. We now define a multiscale decomposition for the in¬ 
tegration in the second regime. Proceeding as we did in the first regime, we define 
1)2 := |'log 2 (K 2 e^)], for h G {f} 2 , • • • , f}i} and h' G { 1)2 + 1, • • • , f}i}, we define 

Uh,uj{^) ■= Xo(2“^||k - p'J^^olln), := Uh',w{^) - /^h'-i,<^(k) 

^ suppf^h,uj, := 6/3,L H SUppf^h',ui- 

The analogues of (4.12), and (4.13) hold with 

ti to 

Ah-iM := i(k) + /^;,_i,^(k)wf-^)(k) + ^ TPi'^'^(k) + wPik). (4.18) 

h'=h h'=l)i 

5 - Intermediate regime. The intermediate region between regimes 2 and 3 
is defined in analogy with that between regimes 1 and 2: we let 

/i,2,to(k) := xo(2“'’"||k-p)(_o||n) - Y1 Xo(2“^"||k - p‘){_Jiii) 

je{o,i,2,3} (4.19) 

/=st2,(.o,i(k) := Xo(2“'’^||k{,j||iii) 

where i )2 '■= [log 2 (K 2 e^)J from which we define gt, 2 ,ui{k) and 5 ,g(j 2 ,a;(k) in the same way 
as in (4.13) with 

ti to 

4,Jk) := i(k) + 4"'4) + wPik). (4.20) 

h'=t)2 h'=l)i 

The analogue of (4.12) holds here as well. 

6 - Third regime. There is an extra subtlety in the third regime: we will see 
in section 9 that the singularities of the dressed propagator are slightly different from 
those of the bare (i.e. non-interacting) propagator: at scale h the effective Fermi points 
p^j with j = 1, 2, 3 are moved to p^j^\ with 

“ PF,illiii ^ (const.) |17|e^. (4.21) 

The central Fermi points, j = 0, are left invariant by the interaction. For notational 
uniformity we set p^q^^ = p‘po- Keeping this in mind, we then proceed in a way 
reminiscent of the first and second regimes: let 1)^ := \}og 2 {'^/13 )\, for /i G {f)/ 3 , • • • , ^ 2 } 
and h' G {f )/3 -|- 1, • • • , \) 2 }, we define 

Uh,Lo,j{k) := xo(2"''||k - lliii)> fh',cj,j(k) := f^h',Lo,jik) - f^h'-i,Lj,j(k) 

^ suppf^h,cu,j, := 8/3,1 n suppf^h',i^,j 

(4.22) 
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and the analogues of (4.12), and (4.13) hold with 

^-i,^,j(k) := i(k) + /^;,_i,,,,(k)#f-')(k) 

h i)i So 

+li-f'(k) + 5: .if''(k)+ 5: lif'tk), 

h'=h h'=ti2 h'=hi 

(4.23) 

7 - Last scale. Recalling that |A:o| ^ vr//3, the smallest possible scale is f )/3 := 
[log2(7r//3)J. The last integration is therefore that on scale h = ()/s + 1, after which, we 
are left with 

< (4.24) 

The subsequent sections are dedicated to the proof of the fact that both F(,^ and 
are analytic in U, uniformly in L, /3 and e. We will do this by studying each regime, one 
at a time, performing a tree expansion in each of them in order to bound the terms of 
the series (see section 5 and following). 


4.2 x-space representation of the effective potentials 

We will compute the truncated expectations arising in (4.4), (4.5), (4.8) and (4.9) 
using a determinant expansion (see lemma 2.1) which, as was mentioned above, is only 
useful if the propagator and effective potential are expressed in x-space. We will discuss 
their definition in this section. We restrict our attention to the effective potentials 
since, in order to compute the two-point Schwinger function in the regimes we are 
interested in, we will not need to express the kernels of in x-space. 

1 - Ultraviolet regime. We first discuss the ultraviolet regime, which differs 
from the others in that the propagator does not depend on the index oj. We write 
in terms of its kernels (anti-symmetric in the exchange of their indices), defined as 

00 - 

V''‘)(V.<«>)=:E 7 »iAn 5 w E E 

'-1 o=(Ql,-,a 2 i) 

ki—k2H-hk2;_l—k2!=0 


'Cki,ai 'Ck2,a2 


The x-space expression for is defined 


as 


V’: 




■ /?|A| 




k,Q; 


^k2;_l,02!-l ^k2!,a2I ' 

(4.25) 


(4.26) 


kee 


i^h) 
13,L 
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so that the propagator’s formulation in x-space is 


5h(x-y):=^ y^ghik) 


(4.27) 


keB 


i^h) 

'f3,L 


and similarly for g^h, and the effective potential (4.25) becomes 

oo „ 

VW(V;(<^)) = EE J dxi • • • dx2l - X2Z, • • • , X2Z_1 - X2/)- 


^Xi,Q;i V^X2,a2 rX2i_i,Q;2i-l rX2Z,Q^2Z 

(4.28) 


with 




^2Z,a(ui>--- >U2Z-l) : = 


(/?|A|) 


2Z-1 


(4.29) 




Remark: From (4.25), W^i^^Ck.) is not defined for kj 0 however, one can easily 

check that (4.29) holds for any extension of to thanks to the compact support 

properties of in momentum space. In order to get satisfactory bounds on 
that is in order to avoid Gibbs phenomena, we define the extension of 14^2^^ (k) similarly 
to (4.25) by relaxing the condition that is supported on and iterating (4.4). 

In other words, we let IF 2 f^(k) for k G be the kernels of defined inductively 

by 

oo 

-/ 3 |A|ez.-V*W('k) := J] + xj,); AT) (4.30) 

Af=l 

in which {'I'k,a}ke 6;3 i,«e .4 is a collection of external fields (in reference to the fact that, 
contrary to they have a non-compact support in momentum space). The use of 

this specific extension can be justified ab-initio by re-defining the cutoff function y in 
such a way that its support is M, e.g. using exponential tails that depend on a parameter 
in such a way that the support tends to be compact as goes to 0. Following this 
logic, we could first define IF using the non-compactly supported cutoff function and then 
take the —)■ 0 limit, thus recovering (4.30). Such an approach is dicussed in [BM02]. 

From now on, with some abuse of notation, we shall identify with and denote 
them by the same symbol which is justified by the fact that their kernels are (or 
can be chosen, from what said above, to be) the same. 

2 - First and second regimes. We now discuss the first and second regimes 
(the third regime is very slightly different in that the index u is complemented by an 
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extra index j and the Fermi points are shifted). Similarly to (4.25), we define the kernels 
of V: 


VW(^(<^)) =: 

1=2 


m\) 


21 - 


tE e 




(ki,- 

ki-k 2 + 


■■+k 2 i_i—k 2 j =0 


, ^2l-l)- 


^ki,Qi,a;i ^k 2 ,0125^2 ^k' 


2i-l5<^2Z-l5^2Z-l 




Note that the kernel 14^2^^ independent of 


i^h)- 
^21,0^21,^^21 ' 

(4.31) 


where x • • • x ■ 

w, which can be easily proved using the symmetry uii i—)• — w*. The x-space expression 


for if)- 


(^h)± . 

k,Q:,a; 


IS 


V^x,Q,a; 


)9|A| 


E 


±i(k-p“ (,).x 7(^/i)± 
Xk,a,(.j • 


(4.32) 


keB' 


0,L 


Remark: Unlike the ^ index in ip^oA i® redundant. Keeping track of 

this dependence is required to prove properties of FF 2 z(k) and 5 /i(k) close to p'^g while 
working in x-space. Such considerations were first discussed in [BG90] in which ^px,a,Lo 
were called quasi-particle fields. 


We then define the propagator in x-space: 
5h,c^(x-y) := ^ Y 


(4.33) 


keB 


P,L 


and similarly for g^h,uj, and the effective potential (4.31) becomes 

OO 

vW(^(<^)) = EEm/ dxi • • • dx2i 1k4^^^^(xi - X2Z, • • • , X2Z_1 - X2z)- 
1=2 ^,o 


and 


in which 


</^X2,«2,‘<J2 rx2i_i,a2l-l ^^ 21 , 021,^^21 

(4.34) 

a;,(Q:,a^) 


,u 2 Z-i) 


T,E|Li(-ihP^o 
(/3|A|)2^-i 


^ ed^f='(-^)Hk2-p7,oy'^i)w^P^{ku • • • , k2z_i). 


(ki,-,k2,_i)eB^(-i 


(4.36) 
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As in the ultraviolet, the definition of is extended to by defining it as the 

kernel of 

-/3|A|e;,-V*WW := ^ + ^);iV) (4.37) 

Af=l 

in which {^'k,a}ke 6/3 i,oe .4 is a collection of external fields. The definition (4.36) suggests 
a definition for (see (4.14) and (4.18)): 

A,Ux):=^ (4.38) 

' keB;3,i 


3 - Third regime. We now turn our attention to the third regime. As discussed 
in section 4.1, in addition to there being an extra index j, the Fermi points are also shifted 
in the third regime. The kernels of V and Q are defined as in (4.31), but with u replaced 
by The x-space representation of dehned as 

{^h)± __ 1 

• olAl Z_^ ’ vk,Q,ajJ 


^S!Z-=m S (4.39) 

and the x-space expression of the propagator and the kernels of V and Q are defined by 
analogy with the first regime: 




(4.40) 




and 


- ,U 2 «-i) :=-- 


^ ,k2,_r). 

{■ku- M2l-l)£l3f l^ 


In addition 


' keB^.L 


(4.41) 

(4.42) 
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4.3 Estimates of the free propagator 

Before moving along with the tree expansion, we first compute a bound on in 
the different regimes, which will be used in the following. 

1 - Ultraviolet regime. We first study the ultraviolet regime, i.e. h G 
{I,... ,M}. 

1-1 - Fourier space bounds. We have 

i(k)-:=-(.W + H„W)-=(l + 

and 

l5h(k)| = |/ft(k)i“^(k)| ^ (const.) 2~^, 
where | • | is the operator norm. Therefore 

l5ft(k)| ^ (const.). (4.43) 

Furthermore, for all mo + ruk ^ 7 (we choose the constant 7 in order to get adequate 
bounds on the real-space decay of the free propagator, good enough for performing the 
localization and renormalization procedure described below; any other larger constant 
would yield identical results), 

|2'^"°C^r5h(k)| ^ (const.) 2-^ (4.44) 

in which denotes the discrete derivative with respect to ko and, with a slightly abusive 
notation, dk the discrete derivative with respect to either ki or k 2 ■ Indeed the derivatives 
over k land on ikQA~^, which does not change the previous estimate, and the derivatives 
over fco either land on fh, l/{iko), or ikoA~^, which yields an extra 2~^ in the estimate. 

Remark: The previous argument implicitly uses the Leibnitz rule, which must be used 
carefully since the derivatives are discrete. However, since the estimate is purely di¬ 
mensional, we can replace the discrete discrete derivative with a continuous one without 
changing the order of magnitude of the resulting bound. 


1-2 - Configuration space bounds. We now prove that the inverse Fourier 
transform of 


5h(x) := 


1 


/3|A| 


E 


^ —2kx 


4(k) 


(4.45) 


kez? 


13 ,L 


satisfies the following estimate; for all mo -|- ruk ^ 3, 

[ dx x^^x^A9h{-^)\ ^ (const.) 


(4.46) 
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where we recall that f dx is a shorthand for dtJ^xeA- Indeed, note that the right 
side of (4.45) can be thought of as the Riemann sum approximation of 




g-ik-X 


5ft (k) 


(4.47) 


where Aqo = {tiGi+t 2 G 2 : ti G [0,1)} is the limit as L —)■ oo of A, see (2.4) and following 
lines. The dimensional estimates one finds using this continuum approximation are the 
same as those using (4.45) therefore, integrating (4.47) 7 times by parts and using (4.44) 
we find 


l5ft(x)| ^ 


(const.) 

1 + (2'*|xo| + \x\y 


so that by changing variables in the integral over xq to 2^xq, and using 


r 

J ^^TTJu[TW 


< (const.) 


we find (4.46). 


J.' 11 1 tigiiiit;. 


2-1 - Fourier space bounds. From (3.8) we 


find 


| 5 ft,a;(k)| ^ (const. 


-t—h 


therefore 

i \9h,M\ ^ (const.) 2'^^ (4.48) 

and for m ^ 7, 

^ (const.) 2-^ (4.49) 

in which we again used the slightly abusive notation of writing to mean any derivative 
with respect to ko, ki or k 2 - Equation (4.49) then follows from similar considerations as 
those in the ultraviolet regime. 

2-2 - Configuration space bonnds. We estimate the real-space counterpart 

of 5ft,a;; 

5 ft,a;(x) ;= ^ 
and find that for m ^ 3, 


j dx \x^gh,i^{x)\ ^ (const.) 2 

which follows from very similar considerations as the ultraviolet estimate. 


(4.50) 
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3 - Second regime. We treat the second regime, i.e. h G {f )2 + 1, • • • , f)i} in 
a very similar way (we skip the intermediate regime which can be treated in the same 
way as either the first or second regimes): 

^ E < (consto 2'“+'“ (4,51) 

and for all mo + mk ^ 7, 

^ (const.) 2-^ (4.52) 

where he := log 2 (e). Therefore for all mo + mk ^ 3, 

J dx ^ (const.) ( 4 , 53 ) 


4 - Third regime. Finally, the third regime, i.e. h G {f )3 + 1, • 
^ l 5 h,i.(k)| ^ (const.) 


^2}: 


/3|A| 


keB 


P,L 


and for all mo + mk ^ 7, 

^ (const.) 2 -^ 

Therefore for all mo + mk ^ 3, 

[ dx \x^°x^'^gh,u^j{^)\ ^ (const.) 


where 


5^/1,oj,j •— 


lAI 


/3|A| 




keB 


13,L 


(4.54) 

(4.55) 

(4.56) 


5 Tree expansion and constructive bounds 

In this section, we shall dehne the Gallavotti-Nicolo tree expansion [GN85], and 
show how it can be used to compute bounds for the Zh, and dehned 

above in (4.4) and (4.8), using the estimates (4.46), (4.50), (4.53) and (4.56). We 
follow [BG90, GMOl, GMIO]. We conclude the section by showing how to compute the 
terms in that are quadratic in Jk^o, from and gh- 

The discussion in this section is meant to be somewhat general, in order to be 
applied to the ultraviolet, first, second and third regimes (except for lemma 5.2 which 
does not apply to the ultraviolet regime). 
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5.1 Gallavotti-Nicolo Tree expansion 

In this section, we will define a tree expansion to re-express equations of the type 

(5.1) 

Af=l 


for h G {/i2) • • • ) “ 1 } (ill the ultraviolet regime /i2 = h\ = M] in the first = (fi, 

h\ = f)o; in the second /12 = i)2, h\ = and in the third, = fl/3, h\ = 1)2), with 


l=q ^ 

='^Yj [ • • • V’S!-0^2i-lV’Sfei 

1 _ r\ ^ 'J 


(5.2) 


Z=0 ro 


(g = 1 in the ultraviolet regime and g = 2 in the first, second and third) in which ro and 
X are shorthands for (tz7i, • • • , 1 x 121 ) and (xi, • • • ,X 2 z); w denotes a collection of indices; 
{a,Lj) in the first and second regimes, {a,Ld,j) in the third, and (a) in the ultraviolet; 
and 11^2^^ (x) is a function that only depends on the differences Xj — Xj. The propagator 
associated with will be denoted fl'(/i+i),(ro,ro')(x — x') and is to be interpreted as the 
dressed propagator g[h+i,uj),{a,a') iii the first and second regimes, and as g{h+i,Lu,j),{a,a') 
in the third. Note in particular that in the first and second regimes the propagator is 
diagonal in the (x indices, and is diagonal in {x,j) in the third. In all cases, we write 

g{h+i),{^,z^'){^-^') = , (5.3) 

' keB;3.L 


where should be interpreted as 0 in the ultraviolet regime, as p'^Q in the first and 

second, and as in the third, see (4.21). 

Remark: The usual way of computing expressions of the form (5.1) is to write the right 
side as a sum over Feynman diagrams. The tree expansion detailed below provides a 
way of identifying the sub-diagrams that scale in the same way (see the remark at the 
end of this section). In the proofs below, there will be no mention of Feynman diagrams, 
since a diagramatic expansion would yield insufficient bounds. 


We will now be a bit rough for a few sentences, in order to carry the main idea 
of the tree expansion across: equation (5.1) is an inductive equation for the which 
we will pictorially think of as the merging of a selection of N potentials via a 

truncated expectation. If we iterate (5.1) all the way to scale then we get a set of 
merges that fit into each other, creating a tree structure. The sum over the choice of 
N's at every step will be expressed as a sum over Gallavotti-Nicolo trees, which we will 
now define precisely. 
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Given a scale h G — and an integer N ^ 1, we define the set 

of Gallavotti-Nicolo (GN) trees as a set of labeled rooted trees with N leaves in the 
following way. 

• We define the set of unlabeled trees inductively: we start with a root, that is 
connected to a node vq that we will call the first node of the tree; every node 
is assigned an ordered set of child nodes, vq must have at least one child, while 
the other nodes may be childless. We denote the parent-child partial ordering by 
v' < V {v' is the parent of v). The nodes that have no children are called leaves or 
endpoints. By convention, the root is not considered to be a node, but we will still 
call it the parent of uq. 

• Each node is assigned a scale label h' G {h + 1, - ■ ■ ,hl + l} and the root is assigned 
the scale label h, in such a way that the children of the root or of a node on scale 
h' are on scale h' + 1 (keep in mind that it is possible for a node to have a single 
child). 

• The leaves whose scale is ^ h\ are called local. The leaves on scale h'l + 1 can 
either be local or irrelevant (see figure 5.1). 

• Every local leaf must be preceeded by a branching node, i.e. a node with at least 
two children. In other words, every local leaf must have at least one sibling. 

• We denote the set of nodes of a tree r by ^(r), the set of nodes that are not leaves 
by TJ(t) and the set of leaves by C:(r). 

Remark: Local leaves are called “local” because those nodes are usually applied a 
localization operation (see e.g. [BG95]). In the present case, such a step is not needed, 
due to the super-renormalizable nature of the first and third regimes. 

Every node of a Gallavotti-Nicolo tree r corresponds to a truncated expectation 
of effective potentials of the form (5.1). If one expands the product of factors of the form 
(V'xfro^^ + V’x^ro^^^) in every term in the right side of (5.1), then one finds a sum over 
choices between and for every (x, tu, ±). We will express this sum as a sum 

over a set of external field labels (corresponding to the labels of which are called 

external because they can be factored out of the truncated expectation) defined in the 
following way. Given an integer Iq q, whose purpose will become clear in lemma 5.2 
(we will choose Iq to be = 1 in the ultraviolet regime, and = 2, 3, 2 in the first, second, 
third infrared regimes, respectively), a tree r G 7)^^^ whose endpoints are denoted by 
(ui, • • • , vn), as well as a collection of integers l.^ := {fii, ■ ■ ■ , Ivj^) G such that 0 q 
and, if Vi is a local leaf, < (.q (in particular, li Iq = q there are no local leaves), we 
introduce an ordered collection of fields, i.e. triplets 

F = ((xi,Ct7i,-|-), {X2,W2, -),••• , (x 2 L_l,ro 2 X,_l, +), {yi2L,^2L, “)) • (5.4) 
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Figure 5.1: example of a tree on scale h up to scale + 1 = h + 5 with 11 leaves, 5 of 
which are local and 6 irrelevant. Local leaves are represented as empty circles, whereas 
irrelevant leaves are represented as full circles. 
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where L + • • • + . We then define the set of external field labels of each endpoint 

Vi as the following ordered collections of integers 

d^vi ■— ( 1 ) ■ ■ ■ ) ‘^Ivi ) ) ■ ■ ■ ) Ivpf ■— {“^Ivi^-i + !,••• , ‘^Ivn') • 

We define the set Vr,i^/Q of external field labels compatible with a tree r G 7^^^ as the 
set of all the collections P = {Pi;}i,e 2 J(r) where Py are themselves collections of integers 
that satisfy the following constraints: 

• For every v G 5J(r) whose children are (r;i, • • • , Vg), Py C Py^ U • • • U Py^ in which, 
by convention, if Vi is an endpoint then Py^ = and the order of the elements of 
Py is that of Pyj through Py^ (in particular the integers coming from Py-^ precede 
those from Py^ and so forth). 

• For all V G 73(r), Py must contain as many even integers as odd ones (even integers 
correspond to fields with a —, and odd ones to a +). 

• If u has more than one child, then Py 7 ^ Py/ for all v' y v 

• For all V G ^(t) \ {uq} which is not a local leaf, the cardinality of Py must satisfy 

|P^| ^ 24 . 

Furthemore, given a node v whose children are (ui, • • • , Vg), we define Ry := Ui=i PvAPv- 

We associate a value to each node v of such a tree in the following way. If u is a 
leaf, then its value is 

(5.5) 

where |Py| denotes the cardinality of Py, and Wy and Xy are the field labels (i.e. elements 
of F) specified by the indices in Py. If u is not a leaf and Ry then its value is 

;= aTyWg{hy),i j dPrAA (5-6) 

T„eT{Ry) ier„ ^ r„eT{Ry) 

where T(Ry), g(hy),h dPTy{i) and are defined as in lemma 2.1 with g replaced by 

ghy, and if the children of v are denoted by (ui, • • • , Vg), then Ry := {PyAPv, ■ ■ ■, PvA^v)- 
If V is not a leaf and Ry = then it has exactly one child and we let its value be /Oy = 1. 
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Remark: The sum over P G 'Pr,i^,eQ is a sum over the assignement of Py for nodes 
that are not endpoints. The sets ly are not summed over, instead they are fixed by 
Furthermore, if Vr,i ,4 = ® (®-§- if ^0 = 9 and r contains local leaves), then the sum 
should be interpreted as 0. 

By injecting (5.6) into (5.7), we can re-write 

-EY. E E C’ n 

N=l T&T{t) DeQJ(T) 

(5.8) 

where T(r) is the set of collections of {Ty G T(R^))„gQj(T-). Moreover, while pi was 
defined in (5.6) if u G ^{t), it stands for if u G C:(r) (note that in this case Ty = 0). 

Idea of the proof: The proof of this lemma can easily be reconstructed from the 
schematic description below. We do not present it in full detail here because its proof 
has already been discussed in several references, among which [BG95, GMOl, GilO]. 


The lemma follows from an induction on h, in which we write the truncated 
expectation in the right side of (5.1) as 


/ dxi • • • dxiv ) • • • 


h,--- ,In 


■£l+i • • • 

In 

n ('^(=£^)+ _|_.^(^+ 1 )+ 

WJV,2J-1'' ^XN,2j,T^ 
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which yields a sum over the choices between and with each choice corre¬ 
sponding to an instance of each “creates” the element {'x.,w,e) in P^. The 

remaining truncated expectation is then computed by applying lemma 2.1. Finally, the 
with Ij < Iq are left as such, and yield a local leaf in the tree expansion, the 
others are expanded using the inductive hypothesis. 

Remark: For readers who are familiar with Feynman diagram expansions, it may be 
worth pointing out that a Gallavotti-Nicolo tree paired up with a set of external field 
labels P represents a class of labeled Feynman diagrams (the labels being the scales 
attached to the lines, or equivalently to the propagators) with similar scaling properties. 
In fact, given a labeled Feynman diagram, one defines a tree and a set of external field 
labels by the following procedure. For every h, we define the clusters on scale h as the 
connected components of the diagram one obtains by removing the lines with a scale 
label that is < h. We assign a node with scale label h to every cluster on scale h. The set 
Pv contains the indices of the legs of the Feynman diagram that exit the corresponding 
cluster. If a cluster on scale h contains a cluster on scale h + 1, then we draw a branch 
between the two corresponding nodes. See hgure 5.2 for an example. 

Local leaves correspond to clusters that have few external legs. They are consid¬ 
ered as “black boxes”: the clusters on larger scales contained inside them are discarded. 

A more detailed discussion of this correspondence can be found in [GMOl, sec¬ 
tion 5.2] among other references. 
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V3 




Figure 5.2: Example of a labeled Feynman diagram and its corresponding tree. Three 
clusters, denoted by vi, V 2 and V 3 , on scale hi, /12 and /14 respectively, are explicitely 
drawn as dotted ellipses. There are 4 more clusters (2 on scale hi, 1 on scale /12 and 1 
on scale / 13 ) which are not represented. The scales are drawn in different colors (color 
online): red for / 15 , orange for / 14 , yellow for / 13 , green for /i 2 and blue for hi. 
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5.2 Power counting lemma 

We will now state and prove the power counting lemma, which is an important 
step in bounding the elements in the tree expansion (5.8) in the first, second and third 
regimes. 

In the following, we will use a slightly abusive notation: given x = (xi,... ,x„), 
we will write x™' to mean “any of the products of the following form” 


^jl,h ■ ■ ■ 

where G {0,1, 2} indexes the components of x and G {1, • • • ,n} indexes the compo¬ 
nents of X. We will also denote the translate of x by y by x — y = (xi — y,..., x„ — y). 
Furthermore, given x™, we dehne the vector m whose i-th component is the number of 
occurrences of x.^i in the product x”^ (note that mo -|- mi -|- m 2 = m). 


The power counting lemma will be stated as an inequality on the so-called beta 
function of the renormalization group flow, defined as 




- IT; 


ifi>« 




(5.9) 


In terms of the tree expansion (5.7), is the sum of the contributions to 11 ^ 2 ^^ whose 
field label assignment P is such that every node v G TJ(r) \ {uq} that is connected to 
the root by a chain of nodes with only one child satisfies |P(,| > 21. We denote the set of 
such field label assignments by 'Pr,z^,£o given r, and In other words, 

contains all the contributions that have at least one propagator on scale h -|- 1. If Z < g, 
then all the contributions have a propagator on scale h -|- 1 , so B 21 = IF 2 Z- 


I Lemma 5.2 

Assume that the propagator ~ '^) can be written as in (5.3). Given h G 

{h* 2 , ■ ■ ■ , — 1}, if Vm G {0,1, 2 , 3}, and 

J dx \x"^gh'{x)\ ^ 

< . , W G {h + 1,--- ,hl}, (5.10) 

where Cg, Ck, Cg and Cq are constants, independent of h, and ^h'ilR) is a shorthand for 

^”^0 2 !™! 2~^'i'^omo+dimi +d2m2) 

in which Ao,Ai,A 2 > 0 , do,di,d 2 ^ 0 , and m* is the number of times any of the xy* 
appears in x™; if 

4 > (5.11) 

Ck Cg 


51 



and 

|(x-x 20’"w^it2u)| ^ 

V/i' G {/i + l,--- ,h\} 

(5.12) 

where (7 ^ ^ for /i' < /i^, Z ^ g for /i' = h\ (in particular, if g ^ £ 0 , then h' = /i^), 
and C 2 Z are constants, then 


1 

M 


dx (x - X 20 ™i?gL(x)| ^ 2'*('^'=-('='=-'^«)')(5;,(m)(C'3C'c^) 'E E E 

iV=l,g7GM 


Pe^r,z^,^0 \i^G2?(r) / \i;Ge(r) / 

iPuQ 1=2/ 

(5.13) 

where Ci, C2 and C3 are constants, independent of Cg, Ck, Cg, Cq and h. 


Remarks: Here are a few comments about this lemma. 

• Combining this lemma with (5.9) yields a bound on fk 2 fL(—)• particular, if 
I ^ io and h < hi, then 

00 

<ix |(x-X2,)’"l4';^(x)|<2'‘l'‘-l'‘-=.»Jft(m)(C3Co‘)' 

^ cf(CgC'G^)^"M n j Y[ {C2CGy^(t2iJur^^^’^^-^^ 

PepT,z.^/o \i’G5?(t) j \'i;G(£(T) 

iPuQ |—2/ 

(5.14) 

• The lemma cannot be used in this form in the ultraviolet regime, since in that case 
the right side of (5.11) is infinite, because Ck = Cg = 1. In the ultraviolet we will 
need to re-organize the tree expansion, in order to derive convergent bounds on 
the series, as discussed in section 6 below. 

• The lemma gives a bound on the m-th derivative of I^ 2 f 2 (—)’ we will need 

in order to write the dominating behavior of the two-point Schwinger function as 
stated in Theorems 1.1, 1.2, 1.3; however, we will never need to take m larger than 
3, which is important because the bound (5.13), if generalized to larger values of 
m, would diverge faster than m! as m —)• 00 . 

• Recall that the propagator appearing in the statement should be interpreted as 
the dressed propagator gh in the first, second and third regimes. Since gn depends 
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on W 21 ^ for h' ^ /i, we will have to apply the lemma inductively, proving at each 
step that the dressed propagator satisfies the bounds (5.10). 

• Similarly, the bounds (5.12) will have to be proved inductively. 

• In this lemma, the purpose of (.q, which up until now may have seemed like an 
arbitrary definition, is made clear. In fact, the condition that (.q > Ckl{ck — Cg) 
implies that Ck — {ck — Cg)\Pv\/2 < 0, Vu G 23(t) \ {uq}. If this were not the case, 
then the weight of each tree r could increase with the size of the tree, making the 
right side of (5.13) divergent. 

• The combination Ck — {ck — Cg)\Py\/2 is called the scaling dimension of the cluster 
V. Under the assumptions of the lemma, the scaling dimension is negative, Vu G 
5J(r) \ {uq}. The clusters with non-negative scaling dimensions are necessarily 
leaves, and condition (5.12) corresponds to the requirement that we can control 
the size of these dangerous clusters. Essentially, what this lemma shows is that 
the only terms that are potentially problematic are those with non-negative scaling 
dimension. This prompts the following definitions: a node with negative scaling 
dimension will be called irrelevant, one with vanishing scaling dimension marginal 
and one with positive scaling dimension relevant. 

• We will show that in the first and third regimes Ck = 3 and Cg = 1, so that the 
scaling dimension is 3 — |P,;|. Therefore, the nodes with |Ey| = 2 are relevant 
whereas all the others are irrelevant. In the second regime, Ck = 2 and = 1, so 
that the scaling dimension is 2 — \P^\/2. Therefore, the nodes with \Py\ = 2 are 
relevant, those with |P^| = 4 are marginal, and all other nodes are irrelevant. 

• The purpose of the factor ^h{lR) is to take into account the dependence of the 
order of magnitude of the different components ko, ki and k 2 in the different 
regimes. In other words, as was shown in (4.46), (4.50), (4.53) and (4.56), the 
effect of multiplying g by xj^i depends on i, which is a fact the lemma must take 
into account. 

• The reason why we have stated this bound in x-space is because of the estimate 
of det(G^^’'’^“^) detailed below, which is very inefficient in k-space. 


Proof : The proof proceeds in five steps: first we estimate the determinant appear¬ 
ing in (5.6) using the Gram-Hadamard inequality; then we perform a change of variables 
in the integral over x^- in the right side of (5.8) in order to re-express it as an integration 
on differences x* — xj; we then decompose (x — X 2 i)™'; and then compute a bound, which 
we re-arrange; and finally we use a bound on the number of spanning trees T(r) to 
conclude the proof. 

1 - Gram bound. We first estimate | det 
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1-1 - Gram-Hadamard inequality. We shall make use of the Gram- 
Hadamard inequality, which states that the determinant of a matrix M whose compo¬ 
nents are given by Mij = Aj © Bj where (Aj) and (Bj) are vectors in some Hilbert 
space with scalar product 0 (writing M as a scalar product is called writing it in Gram 
form) can be bounded by 

I det(M)| ^ \/Aj 0 Aji/Bj 0 B*. (5.15) 

i 

The proof of this inequality is based on applying a Gram-Schmidt process to turn (Aj) 
and (Bj) into orthonormal families, at which point the inequality follows trivially. We 
recall that is an (n^, — (s„ — 1 )) x (n^, — (s„ — 1 )) matrix in which Sy denotes 

the number of children of v and if we denote the children of u by (ui, • • • , ), then 

riv = \Rv\/2 = {'^i=i\Pvi\ — |T’i,|)/2. Its components are of the form te9(hy),e (see 
lemma 2.1), with t(j= Uj • uj in which the ttj are unit vectors. 


1-2 - Gram form. We now put ~ Gram form 

by using the k-space representation of gh in (5.3). Let Ti = x b, d, b}) denote 

the Hilbert space of square summable sequences indexed by (k, a) G 13/3,l x {a, 6, a, 6}. 
For every h G {h^, • • • , /i| — 1} and (x, a) G ([0, /S) x A) x {a, b, d, b}, we define a pair of 
vectors (Ai^^(x), Ba^^(x)) G by 




(5.16) 


where A^^^(k) denotes the a-th eigenvalue of 5 j)(k) 5 /i(k) (i.e. the singular values of 
ghO^)) and V^P(k) and UP\'k) are unitary matrices that are such that 

^;,(k) = l/('*)t(k)z)W(k)t/('^)(k), 

where .D(^)(k) is the diagonal matrix with entries Ai^^(k). We can now write gu as 

5(h),(a,a') (x - x) = A^'")(x) © B^ (x ) (5.17) 

where © denotes the scalar product on %. Furthermore, recalling that | 5 fe(k)| is the 
operator norm of 5 /i(k), so that |5h(k)| = max spec-^ gj^(k)gh(k), we have 

Ai'“l(x)0Ai'“l(x) = Bi'‘l(x)OBi'“l(x)<^ |9,.(k)| < Cg2‘'=‘-'«>'‘ (5.18) 




The Gram form for is then 

- Xj) = {Ui ■ Uj)(A^,(Xj) © Byy.{Xj)) 


(5.19) 
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(5.20) 


so that, using (5.15) and (5.18), 

I det I ^ 

2 - Change of variables. We change variables in the integration over x.,.. For 

every v G ^^(t), let Pv =: (i|^\ • • • We recall that a spanning tree T G T(t) is 

a diagram connecting the fields specified by the I^’s for v G 2:(t): more precisely, if we 
draw a vertex for each v G ^(t) with |/„| half-lines attached to it that are labeled by the 
elements of then T G T(r) is a pairing of some of the half-lines that results in a tree 
called a spanning tree (not to be confused with a Gallavotti-Nicolo tree) (for an example, 
see figure 5.3). The vertex Vr of a spanning tree that contains the last external field, i.e. 
that is such that G is defined as its root, which allows us to unambiguously 
define a parent-child partial order, so that we can dress each branch with an arrow that 
is directed away from the root. For every v G ^{r) that is not the root of T, we define 

G as the index of the field in which T enters, i.e. the index of the half-line of T 
with an arrow pointing towards v. We also dehne ;= Now, for every v G ^(t), 
we define 

“ ^j{«) 

for all G Iv \ and given a line of T connecting to we define 

We have thus dehned (Xlveefr) ^ variables z, so that we are left with Xj{„^), which 

we call xq. It follows directly from the definitions that the change of variables from x.,. 
to {xo, where Ir = UveC(r) has Jacobian equal to 1. 

3 - Decomposing (x —X 2 /)™'. We now decompose the (x —X 2 z)”^ factor in (5.13) 

in the following way (note that in terms of the indices in X 2 l = ^jm)- (x “ ^ 2 l)^ 
is a product of terms of the form {xj^i — Xj(vr) j) which we rewrite as a sum of zj/fis 
for V G ^(t) on the path from to j, a concept we will now make more precise, j 
and are in and respectively, where v{j) is the unique node in C:(r) such 

that j G Iv{j)- There exists a unique sequence of lines of T that links Vr to v{j), which 
we denote by ((ji, j(), • • • , ijp,jp)), the convention being that the line is oriented 

from j to f. The path from to j is the sequence zj^, z^/ ,Zj^, - ■ ■ and so forth, until 
j is reached. We can therefore write 

p 

^j,i ~ 

p=l 

4 - Bound in terms of number of spanning trees. Let us now turn to the 
object of interest, namely the left side of (5.13). It follows from (5.2) and (5.8) that 

<U)=i: E E e n (5.21) 
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Figure 5.3: example of a spanning tree with = 5 and \Pvi \ = \Pv 2 \ = {Pvil = \Pv 4 \ = 4, 
I Fug I = 6; whose root is V 2 - 


Therefore, using the bound (5.20), the change of variables defined above and the decom¬ 
position of (x — X 2 z)™ described above, we find 


1 


^ / dx (x - X 2 /)”^F^!^I^(x 


/?|A| 


E E EE/^^o E 

' ^=i,er^'*)TeT(r) iu F^Pr,l^,io 

\Pvq\=21 


E 

J2i'nie+mv)=m 


E»(r) E"'' 


(^ 1 ) l) 


Yl ( f dze 

ter„ 


t>Sl£(T) 

(5.22) 

(we recall that by definition, if u G ^(t), ly = Py and \Iy\ = 2 ly) in which we inject (5.10) 
and (5.12) to find 


/?|A| 


dx 


(x - X 2 z)”^Fg^(x) ^ EE E (^i^h{rn) 


^=lrer^'*)^6T(r) PS^u.^u/o 
IAjq \ =2l 

_ /~^^V G1 1 ^yhy ((c/g Cq^Tly C^(^Sy l)) 

g I'-'G W ^ 

t'S2J(T) ^ 

(5.23) 
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in which Ci is an upper bound on the number of terms in the sum over {mi) and (m„) 
in the previous equation, and C 2 denotes the number of elements in the sum over 
Recalling that n„ = \Rv \/2 = I-PdJ — \Pv\)/ 2 , we re-arrange (5.23) by using 


[ Y. h,\R,\ = -h\P,,\- Y 1 ^- 1 + 

Y. K{sv-'^) = -h- Yj 1 + Y. ~ 

( ^!S^3(T) vG*XJ(t) VGI£(t) 


and 


|R„| = |I,J - |P, 
^ (s,-l)=iV-l 


Vo I 


to find 
1 


/3|A| 


dx 


x-x20”^<l(x) E EE Y 

N=i tgt(t) R PePr.i,,eg 

IPvqI =21 

.2Hck-(ck-Cs)l)^^(^mj J_2Cfe-(cfe-C3)^ 


^'S2J{T) 


vGS(t) 


(5.24) 


5 - Bound on the number of spanning trees. Finally, the number of choices 
for T can be bounded (see [GMOl, lemma A.5]) 

E n ^ 3 ^"- ( 5 - 25 ) 

reT(r) t>e»(T) 

so that by injecting (5.25) into (5.24), we find (5.13), with C 2 = and C 3 = □ 


5.3 Schwinger function from the effective potential 

In this section we show how to compute in a similarly general setting as 

above; consider 

-m^h - 

= Y {Y^^^\Y^YY^^\Ya);N) 

N=1 

(5.26) 

for h £ {/i- 2 , • • • — This discussion will not be used in the ultraviolet regime, so we 

can safely discard the cases in which the propagator is not renormalized. Unlike (5.1), it 
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is necessary to separate the a indices from the indices, so we write the propagator 

of as g(h+i,-uu),(a,a') where w stands for oj in the first and second regimes, and (cu, j) 
in the third. 

We now rewrite the terms in the right side of (5.26) in terms of the effective 
potential Let 

:= 4a)- (5-27) 

Note that the terms in are either linear or quadratic in Jk,a, simply because the two 
J variables we have at our disposal, 4«i ’ ^^ct 2 ’ Grassmann variables and square to 
zero. We define the functional derivative of with respect to 

:= J #4 vW(V;). 


Lemma 5.3 


Assume that, for h = h\, 

xV'H’p, Aa) = + EfW.C'.t'(kjAo.+ 

cx.' 

a' 


(5.28) 

for some SQ^42(k)) Then (5.28) holds for h G {h 2 ,---,hi — 1} as 

well, with 


:= + E 4L4'4)5(h+i,^),(a»,aO(k) 

,zu 

GE")(k) := G;E'^(k) + E 

a'^ ,zu 

<1;’(k) := dlL'A) - E<571“'A)'LL..„.|(k) 

cl" 

CS’(k) := CA'’(k) - E »l'|L,„.)(k)G;i!2.(k) 


(5.29) 


(5.30) 


and 


ol' ,a" ,uj 


AiAjk) ;=4';S(k)+ E '!LtA’(k)9(ft+i,„),(„',„")(k)v!*„y'(k) 

- E '5:;S'(k)'L!|L.,„")(k)G;,™ (k) 

a' ^a" 

(5.31) 
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in which the sums over a are sums over the indices of g. 


The (inductive) proof of lemma 5.3 is straightforward, although it requires some 
bookkeeping, and is left to the reader. 

Remark: It follows from (4.24) and (2.24) that the two-point Schwinger function S 2 (k) 
is given by S 2 (k) = s*^*’'^^(k) (indeed, once all of the fields have been integrated, = 

J+sW(k)J-). Therefore (5.31) is an inductive formula for the two-point Schwinger 
function. 


6 Ultraviolet integration 

We now detail the integration over the ultraviolet regime. We start from the tree 
expansion in the general form discussed in section 5, with q = Iq = I and h\ = M; 
note that by construction these trees have no local leaves. As mentioned in the first 
remark after lemma 5.2, we cannot apply that lemma to prove convergence of the tree 
expansion: however, as we shall see in a moment, a simple re-organization of it will allow 
to derive uniformly convergent bounds. We recall the estimates (4.46) and (4.43) of gu 
in the ultraviolet regime: for niQ + ruk ^ 3 

J d-K x'^<^x^^\gh{^)\ ^ (const.) 

< 1 _ ( 6 . 1 ) 
2^ lfffe(k)l ^ (const.). 

Equation (6.1) has the same form as (5.10), with 

Cg = Ck = 1 , ^h{m) = 


We now move on to the power counting estimate. The first remark to be made 
is that the values of the leaves have a much better dimensional estimate than the one 
assumed in lemma 5.2. In fact, the value of any leaf, called (x), is the antisymmetric 
part of 

,Ct2 ^^3 (5(xi — X 2 )( 1 (X 3 — X4)l/Wai,a3(xi — X 3 ) (6.2) 

SO that 

^ y dx |(x - X4nri"'(x)| < dilt/l. (6.3) 
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1 - Resumming trivial branches. Next, we re-sum the branches of Gallavotti- 
Nicolo trees that are only followed by a single endpoint: the naive dimensional bound 
on the value of these branches tends to diverge logarithmically as M —)■ oo, but one can 
easily exhibit a cancellation that improves their estimate, as explained below. Consider 
a tree r made of a single branch, with a root on scale h and a single leaf on scale M -\-l 
with value The 4-field kernel associated with such a tree is := 

The 2-field kernel associated with r, once summed over the choices of and over the 
held labels it indexes for /i -|- 1 < /i„ ^ M, keeping and its held labels hxed, can be 
computed explicitly: 



X 


/(5(x 


E 

0-2 ' 


dy 



(6.4) 

If one were to bound the right side of (6.4) term by term in the sum over h' using 
the dimensional estimates on the propagator (see (4.46) and following), one would hnd 
a logarithmic divergence for f q,/)(x)I) i-®- a bound proportional to M — h. 

However, the right side of (6.4) depends on propagators evaluated at xq = 0 (because 
rc(x) is proportional to (5(xo)), so we can use an improved bound on the propagator gh': 
the dominant terms in 5 /i(k) are odd in /cq, so they cancel when considering 


dhOi-)- 


From this idea, we compute an improved bound for |^/i(x)| with xq = 0: 


\gh{0,xi,X 2 )\ ^ 

ki,k2 


'^9h(k) 

ko 


^ (const.) 2 


All in all, we hnd 

j dx |x™A:^|J|^_^,)(x)| ^ ^ J dx \ {x- X4 )^k[^1{x)\ ^ ^^\U\ (6.5) 

for some constant € 4 . We then re-organize the right side of (5.7) by: 

1 . summing over the set of contracted trees which is dehned like but for 

the fact that every node v >~ vq that is not an endpoint must have at least two 
endpoints following it, and the endpoints can be on any scale in [h -|- 2 , M -|- 1 ]; 

2 . re-dehning the value of the endpoints to be with 1 ^ = 1 , 2 . 

2 - Contracted tree expansion. We can now estimate the “contracted tree” 
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expansion, by repeating the steps of the proof of lemma 5.2, thus finding 
1 ^ 


/3|A| 


dx 


(x-x^rw-Wfe) 




E E E E E 

I Pun I 


Cl 


e UQi 

n r 


— 2 


hyi^Su 1 ) 


n 




dS£(t) 


( 6 . 6 ) 

for two constants ci and C 2 in which the sum over is a sum over the A G {1,2}. It 
then follows from the following equation 

Kis, - 1) = hiN - 1) + Y1 

dSOJIt) vS2)(r) 

in which Ny denotes the number of endpoints following v G t, which can be proved by 
induction, that 


1 


/3|A| 


dx 


(x-x2 ,rM4'l(x) 


< 5;(ic/|c3)^2-'‘|«-‘| e E E n 2-1"--'). 

|-f\)Q 1 = 2 / 

(6.7) 

Furthermore, we notice that by the definition of ,, |P„| ^ 2Ny + 2. In particular, 
for V = vq, 21 ^ 2N + 2, so the sum over N actually starts at maxjl, / — I}: 

^^dx |(x-X20™IFi;'^i(x) 

OO 

< E (|c/|c3)^2-'‘|«-‘| e E E n 2-i""-‘i. 

W=max{l,/—1} Lr ^ 

\Pvq\=21 

( 6 . 8 ) 


3 - Bound on the contribution at fixed N. We temporarily restrict to the 
case N > 1. We bound 

E E E n 

A PePr.i.^,! i;e»(T) 

Since Ny ^ 2 and |iA| ^ 2Ny + 2, V// G (0,1), 

-{Ny - 1) ^ min {2 - -l} ^ (1 - ^)min{2 - -l} - ^ ^ -(1 - a)^ - d 
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so that 


E E E n 2 

-r i;e!»(r) 


-(1-/2)^9-M 


3-1 - Bound on the field label assignments. We bound 

PePi-.i^,! DeQJ(r) 

We proceed by induction: if vq denotes the first node of r (i.e. the node immediately 
following the root), (ui, • • • , Vg) its children, and (ri, • • • , Tg) the sub-trees with first node 
(ui, • • • ,Vg), then 


l-Puj |H-l-l-P-Us I 

E n E E E 

PePr.i^.i ^eQJ(r) PieP(Ti) P.eP(r,) pvo=o 


Pvi I + • • • + l-fv 


• 2 2 

i=i ^!e®(Ti) 




=ni E (i+2-V)i^.i 

i=l \PieP(ri) 

SO that by iterating this step down to the leaves, we find 


II 2-(^-^)^ 

vG^in) 


M-h 

\Pv\ ^ I l-fj ■ 


^ H 2-(i-^)^ ^ ^ Cl 

PSPt,!-,! t'SQJfr) \ P=0 / 


for some constant Cp. 


3-2 - Bound on trees. Finally, we bound 

E n 

T-(zrW vg^(t) 

We can re-express the sum over r as a sum over trees with no scale labels that are such 
that each node that is not a leaf has at least two children, and a sum over scale labels: 

E = E E 

r*GT^hGHh{T*) 

in which denotes the set of unlabeled rooted trees with N endpoints and H/i(r*) 
denotes the set of scale labels compatible with t* . Therefore 

^ ^ 2~^ = ^ ^ ^ ^ ‘ 2 ~ 

DSQJfr) heHh(r*) ve»(r*) 
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in which p(v) denotes the parent of v, so that 


n 2 -'- 

T(,r^)v&0{T) T*(iT^v&0{T*)q=l r*eT^ 

for some constant Ct,i, in which we used the fact that |23(r*)| ^ N. Furthermore, it is 
a well known fact that Ylr* ^ ^ [GMOl, lemma A.l], the proof is based on 

constructing an injective map to the set of random walks with 2N steps: given a tree, 
consider a walker that starts at the root, and then travels over branches towards the 
right until it reaches a leaf, and then travels left until it can go right again on a different 
branch). Therefore 

2 “^ ^ ( 6 . 10 ) 

2^e2J(r) 

for some constant Ct- 


< E n E 


c 


N 

T,1 



3-3 - Conclusion of the proof. Therefore, by combining (6.9) and (6.10) 
with the trivial estimate 1 ^ 2^, we find 

-T 


Tat ^ (const.)'^. 


( 6 . 11 ) 


Equation (6.11) trivially holds for N = 
y |(x - 


1 as well. If we inject (6.11) into (6.8) we get; 

OO 

^ Y (|t/|C")^2“'*(^"^) (6.12) 

7V=max{l,(—1} 


for some constant C" and h ^ 0. In conclusion, if |?7| is small enough (uniformly in h 
and 1 ), 


1 


dx 


,X - X2i ) 






for some constant Co > 0. 


(6.13) 


7 First regime 

We now study the first regime. We consider the tree expansion in the general form 
discussed in section 5, with = f)o and q = io = 2, so that there are no local leaves, 
i.e., all leaves are irrelevant, on scale f)o + 1- Recall that the truncated expectation 
in the right side of (5.1) is with respect to the dressed propagator Qh+i in (4.13), so that 
(5.1) is to be interpreted as (4.8). A non trivial aspect of the analysis is that we do not 
have a priori bounds on the dressed propagator, but just on the “bare” one gh,Lo, see 
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(4.50), (4.48). The goal is to show inductively on h that the same qualitative bounds 
are valid for gh,Lj, namely 


j dx |x’”®„{x)| < 

_ I keB,9,r 

which in terms of the hypotheses of lemma 5.2 means 

Cfc = 3, Cg = 1, :Sh{m) = 2-^^. 
Note that Iq = \ck/{ck — Cg)~\ > Ckl{ck — Cg), as desired. 


(7.1) 


7.1 Power counting in the first regime 

It follows from lemma 5.2 and (6.13) that 


1 


/3|A| 


dx 


!x-X2,rBg^,„(x) 


^ 2h{Z—2l) 2—mh 'y 'y y y (j'l^ | ^|max(l,;„ —1) 

oW i;S2J(r) dSi£(t) 


7V=1, 


-(h) I 


-r PeP, 


I PvQ \— 2 l 


for two constants C[ and C". 


(7.2) 


1 - Bounding the sum on trees. First, we notice that the sum over can be 
written as a sum over h, - ■ ■ ,In, so that it can be moved before focus on the 

sum 

E E n 2(3-l^^l). (7.3) 

pp'p(^) 

\Pvq\=21 

We first consider the case I ^ 2. For all 6 G (0,1), 

^ ^ n 


rsT^'"^ 2 i)eQJ(r) 

\Pvq\=21 


\ Pvq \=‘21 


and since io = 2, |Py| ^ 4 for every node v that is not the first node or a leaf, so that 
3 — iPy] ^ —|Py|/4. Now, if N ^ 2, then given r, let v* be the node with at least two 
children that is closest to the root, and h* its scale. Using the fact that |Pi,| ^ 2Z + 2 for 
all V < V* and the fact that r has at least two branches on scales ^ /i*, we have 


2®(3-|Pd) ^ 

dSOJIt) 


2 eht 
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If = 1, we let h* := 0, and note that the same estimate holds. Therefore 


iPuQ |= 2 Z 






20(21 


-l)ih-h*)+ 2 eh* 


h*=h+l 


E E n ^ 

PSPt ,!.^,2 fe 27 ( T )\{ iJo } 


,-( 1 - 0 )^ 


which we bound in the same way as in the proof of (6.11), i.e. splitting 

(1 - 0)1^ = (1 - e)(i -p,)l^ + (1 -> (1 - «)(1 - ,.)L^ + (1 - 9),, 
for all fj, G (0,1) and bounding 

t'eQ7(T)\{i)o} 


and 

^ C^. 

rST^'"^ i'eQJ(T)\{t;o} 

Therefore if ^ ^ 2, then 


E E 


N 


- 6 Tv 


(h) 


PeP, 


(h) 

T,lr,2 


tie®(T)\{iJo} 


p- 


2 = 1 


\Pvq\—21 


(7.4) 


Consider now the case with 1 = 1. If = 1 then the sum over r is trivial, i.e., 1~Y 
consists of a single element, and the sum over P can be bounded as 


E 


2(3 i-p«i) ^ 2^ 


-paV^^^ pG^(r) 
\Pvq \=2 


Per. 


(h) 

t,Zi,2 


\ Pv ()\—2 


PG 51 (r): 

v^v' 


(7.5) 


where v' is, if it exists, the leftmost node such that |Py| > 4, in which case 4 — |P„| ^ 
— |Pu|/3; otherwise, we interpret the product over v as 1. Proceeding as in the case I ^ 2, 
we bound the right side of (7.5) by 


0 

2 hch ^ 2'*C"C''E (7.6) 

h^/ =/i +2 
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If ^ 2, then we denote by t* the subtree with v* : v* as first node, and r' the linear 
tree with root on scale h and the endpoint on scale h* , so that rr' U t* . We split (7.3) 
as 


Af+l -2 


E E E ( n 23-1^’^!) ( E E n (7-7) 

|P„ol=2 \P^*\=2l* 


l*=2 h*=h PeP^,_,*_ 2 : 


The sum in the last parentheses can be bounded as in the case 1^2, yielding h 22 Sh* 
The remaining sum can be bounded as in (7.5)-(7.6) so that, in conclusion, 


-2 


^ ^ ^ ^ 2^^ Tw|) ^ ^ ^ 2^ ^ ^ i229{h' h*)t22dh* 

h*=h h'=h-\-2 

\Pvq\=2 


q~{h) pp-p(^) 

rt !^ x-t/ 


2-1 - ^ = 1. Therefore, if / = 1, (7.2) becomes (we recall that (7 = 2 > 1 so 
that B 2 = IT 2 , see (5.9)) 


/ oo 00 

dx |x”^IT2^tUW| ^ ^ ^ (C"'|t/|)E)Iimax(l,i,-l) (7 9) 

\r_ 17 7 o 


N —1 


Assuming \U\ is small enough and using the subadditivity of the max function, we 
rewrite (7.9) as 

J dx |x™IT2^2,aW| ^ 22 '* 2 -™'^Ci|t/| (7.10) 

which we recall holds for m ^ 3. 

2-2 -1^2. Similarly, if ^ ^ 2, 


/3|A| 


dx 


oo 

^ 2^h(3—2l^29)2^—mh ^ ^ 


E 


'"'Ir/hEili max(l,Zj-l) 


iC'{'\U\) 


N=1 ^lj\i^2 

(/i —1)H-h(^iv —1)^/—H-<5iv,i 

in which the constraint on /i, • • • , /w arises from the fact that, if > 1 , 

1^0 1 ^ |Co|- 2 (A^-l), 


(7.11) 
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while, if = 1, \Pvq\ < \Ivo\- Therefore, assuming that |?7| is small enough and sum¬ 
ming (7.11) over h, we find 


1 


m 

1 


l m 

for Z > 3 and m < 3. 


I <ix |(s-X4r»'W„(x)| 

(x- (x)| < 2''P-2'+2''>2 -’”'>(c,|c/|)‘-‘ 


cZx 


(7.12) 


Remark: The estimates (7.2) and (7.8) imply the convergence of the tree expansion 
(5.8), thus providing a convergent expansion of ^ in U. 


7.2 The dressed propagator 

We now prove the estimate (7.1) on the dressed propagator by induction. We 
recall (4.13) 

&,.4(k))-‘ = 4y(k)i<'*"l(k) (7.13) 

with 

to 

i(^’-)(k) := i(k) + irf'^(k) 

h'=h+l 

whose inverse Fourier transform is denoted by Note that (7.10) on its own does 

not suffice to prove (7.1) because the bound on 


to 

/^4,co(k)lTf)(k)+ ^ Wf)(k) (7.14) 

h'=h+i 


that it would yield is (const.) |Z7| whereas on the support of fh,uj, g~^ ~ 2^, which we 
cannot compare with \U\ unless we impose an e-dependent smallness condition on Z7, 
which we do not want. In addition, even if (7.14) were bounded by (const.) \U\2^, we 
would have to face an extra difficulty to bound g in x-space: indeed, the naive approach 
we have used so far (see e.g. (4.46)) to bound 



would require a bound on d^gh,uj{k) with n > m -|- 3 (we recall that the integral over x 
is 3-dimensional), which would in turn require an estimate on 
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for h' > h, which we do not have (and if we tried to prove it by induction, we would 
immediately find that the estimate would be required to be uniform in n, which we 
cannot expect to be true). 

In order to overcome both of the previously mentioned difficulties, we will expand 

^ (h') 

W 2 ' at first order around contributions up to first order in k — p^g will 

be called the local part of W 2 ■ Through symmetry considerations, we will write the 
local part in terms of constants which we can control, and then use (7.10) to bound the 
remainder. In particular, we will prove that 11^2^^ (p'^g) = 0 from which we will deduce 
an improved bound for (7.14). Furthermore, since the k-dependance of the local part is 
explicit, we will be able to bound all of its derivatives and bound g in x-space. 

1 - Local and irrelevant contributions. We define a localization operator: 

C : Ah,^{-x) I—> (l(x) J dy Ah,uj{y) - d^5{x.) • J dy yAh,uj{y) (7.15) 

where S(x) := d(xo)dxi ,o 42, g and in the second term, as usual, the derivative with respect 
to xi and X 2 is discrete; as well as the corresponding irrelevator: 

7^ :=!-£. (7.16) 

The action of C on functions on k-space is (up to finite size corrections coming from the 
fact that L < 00 that do not change the dimensional estimates computed in this section 
and that we neglect for the sake simplicity) 

TA/i^^(k) = A/i^^(p‘^q) + (k — p^po) • clk^/i,i.j(Pp,o)- ("^•1'7) 

Remark: The reason why C is defined as the first order Taylor expansion, is that its 

" ih') 

role is to separate the relevant and marginal parts of IF 2 ^ from the irrelevant ones. 
Indeed, we recall the definition of the scaling dimension associated to a kernel W 2 (see 
one of the remarks after lemma 5.2) 

Cfc (cfc Cg) = 1 

which, roughly, means that ^ is bounded by = 2 ^'. As was remarked 

above, this bound is insufficient since it does not constrain Ylh'^h ^2 smaller 

than 2^ ~ g~^. Note that, while lF 2 ^^(k) is bounded by 2^', irrespective of k, (k — 
Ppo) ■ ^^ 2 ^ improved dimensional bound, proportional to 2 ^~^'2^', where 

2^ ~ |k — Ppgl; in this sense, we can think of the operator (k — Ppg) • <9k as scaling 
like 2^~^ . Therefore, the remainder of the first order Taylor expansion is bounded 
by ^2^ = 2^^“^ and thereby has a scaling dimension of —1 (with respect to 

h'). Thus, by defining C as the first order Taylor expansion, we take the focus away 
from the remainder, which can be bounded easily because it is irrelevant (i.e., it has 
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negative scaling dimension), and concentrate our attention on the relevant and marginal 
contributions of W 2 . See [BG95, chapter 8] for details. 

We then rewrite (7.13) as 


9h,M = fh,M ' ( 1 +(7^^,^(k)) (Li[;,],^(k))) ' (7.18) 

where Lgj/j] is a shorthand for 

L0[ft],a;(k) := {Uh+1,M - Uh-2,M) (/lih,a;(k)) 

(we can put in the {f^h+i,LoO^) — f!ih- 2 ,uii^)) factor for free because of the initial fh,uiO^))- 
2 - Local part. We first compute CAh^uO^)- 

2-1 - Non-interacting components. As a first step, we write the local 
part of the free inverse propagator as 


where 


£i(k) 


/ iko 

71 

0 

r 

\ 

71 

iko 


0 


0 

r 

iko 

73^ 


V ^ 

0 

73^ 

iko 

/ 




(7.19) 


(7.20) 


2-2 - Interacting components. We now turn to the terms coming from 
the interaction. We first note that ^ satisfies the same symmetries as the initial 
potential V (2.20), listed in section 2.3. Indeed, is a function of V and a quantity 
similar to (2.30) but with an extra cutoff function, which satisfies the symmetries (2.32) 
through (2.38). Therefore 

= wf'^(ii^k) = uiW2^'''^(i?,,k)cJi = -asW^'^'\lk)a3 

=»'f'(Pkf=(; ) <’(r-k) (; 

(7.21) 

^ (h^) 

This imposes a number of restrictions on CW 2 '■ indeed, it follows from propositions F.l 
and F.2 (see appendix F) that, since 

Pf,o = “Pf,o ~ ~ -^/iPf,o = -^Pf,o = 'PPf,o ~ ^Pf,o (7.22) 
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in which R^, R^, I, P and T were defined in section 2.3, we have 


£W^f)(k' + p‘^,o) = - 


/ Kh'ko 0 Uh'^* \ 

lijj-h' iCh'ko Vh'i 0 

0 i^h'C iCh'ko 'jskh'C 

\ Vh'^ 0 'Jskh'C iCh'ko j 


(7.23) 


\Ch'\ ^ (const.) |c/|2^', IC/i'l ^ (const.) |C/|2^', \flh'\ ^ (const.) \U\2'^’^' 
\vh'\ ^ (const.) \U\2'^', \vh'\ ^ (const.) \U\2'^'~'^^. 


with {Ch'■, Jj'h'■, kh'■, Ch'■, h'h') £ Furthermore, it follows from (7.10) that if h' ^ f}o, then 

(7.24) 

St.; |L/ |z, 

Injecting (7.19) and (7.23) into (4.14), we find that 

h 0 VhC* \ 

... t n I 

(7.25) 


where 


Zh ■ = 




/ izhko Jirkh 

0 

VhC* \ 

,a;(k' + Pf,o) — 

- 

'yi'ihh izhko 

0 VhC* 

o 

KKJ' 

* 

0 

lOVhC 



\ VhC 0 

73VhC* 

izhko ) 

to 


to 


to 

1 + Ch', 

rhh 

•= 1 + k'h', 

Vh ■ = 

l+'^kh 

h'=h 


h'=h 


h'=h 


to 


to 


•= 1 + 

E 

, Ch', Vh'-=l + 




(7.26) 


h'=h 


h'=h 


By injecting (7.24) into (7.26), we find 

\mh - 1| ^ (const.) |[/|, \zh - 1| ^ (const.) |t/|, \zh - 1| ^ (const.) \U\, 
\vh - 1| ^ (const.) |17|, \vh - 1| ^ (const.) \U\. 


(7.27) 


2-3 - Dominant part of CAh^u)- Furthermore, we notice that the terms 
proportional to rfih or ai'c sub-dominant; 




(7.28) 


'C^/i,tj(k -|- p^^o) ~ ~ 


(7.29) 


where 

/ izhko 0 0 VhC* \ 

0 izhko VhC 0 

0 VhC* izhko 0 

\ VhC 0 0 izhko J 

Before bounding ai, we compute the inverse of (7.29): using proposition B.l (see ap¬ 
pendix B), we find that if we define 


ko := Zhko, ko := Zhko, C ■= VhC 


(7.30) 
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then 


and 


+ Pf,o) = (koko + |C| 


'^^/i,a;(k)(k^ + Pf,o) — “ 


(fepfco + 1^1^ 
det 


( -iko o_ q r \ 

0 —iko C 0 

0 C -ih 0 

V C 0 0 -iko } 


In particular, this implies that 


l^^/ii(k' + PF,o)l ^ (const.) 2 


which in turn implies 


|iTi(kq| ^ (const.) 2 


he — h 


(7.31) 


(7.32) 


(7.33) 

(7.34) 


3 - Irrelevant part. We now focus on the remainder term TZAh^^{k.) Lg[;j]^^(k) 
in (7.18), which we now show to be small. The estimates are carried out in x space. We 
have 


dx 




= / dx 


dy (Lg[/i],^(x - y) - Lg[;,]_^(x) + y5xLg[;,]_^(x)) 


which, by Taylor’s theorem, yields 


dx 




9 f , 

< - max / dy 

4 *,i J 


ViVj^i^Jiy) 


■ max j I(x) I 


in which we inject (7.10) and (4.49) to find, 

j dx < 2^ (const.) \U\. 

Similarly, we find that for all m ^ 3, 

I dx|x-77wi"J^L0[;,],,(x)| ^2^2—'^(const.) \U\. 

This follows in a straightforward way from 

j dy ynW^'l^\y)hQihl^{x - y) = j dy yVk^tV) (Lg[ft],a;(x - y) - Lg[;,]^^(x)) 


(7.35) 


(7.36) 
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and, for 2 ^ m ^ 3, 

J dy y^nW^^J{y)hgih],u:{^ - y) = ^ dy y'^W^^J{y)hgih],co{^ - y)- 

Remark: The estimate (7.36), as compared to the dimensional estimate without 77, is 
better by a factor \ This is a fairly general argument, and could be repeated with 

1^0[h],a; replaced by the inverse Fourier transform of fh,ui- 


dx 




^ 2 


2h—mh 


(const.) \U\. 


(7.37) 


Finally, using (7.36) and the explicit expression of g, we obtain 


dx 


x"^77i('*’‘^) *L0[;,],^(x) ^2'^2-”^^(const.) (l + |t/||h|). (7.38) 


4 - Conclusion of the proof. The proof of the first of (7.1) is then completed 
by injecting (7.29), (7.34), (7.28), (7.27) and (7.38) into (7.18) and its corresponding 
x-space representation. The second of (7.1) follows from the first. 


7.3 Two-point Schwinger function 

We now compute the dominant part of the two-point Schwinger function for k 
well inside the first regime, i.e. 

f)o-i 

k G := y suppfh,uj- 


Let 

/ik := max{/r : A,^(k) / 0} 
so that if /i 0 {/ik, dk — 1}, then fh,u)0^) = 0. 

1 - Schwinger function in terms of dressed propagators. Since /ik ^ 
the source term +'*^k« 2 '^ko 2 constant with respect to the ultraviolet fields, 

so that the effective source term defined in (5.27) is given, for /i = f)o, by 

Jk,„) = (7.39) 

which implies that is in the form (5.28) with 

gifto) = s('’o)(k) = 0, = 0. 
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Therefore, we can compute for h G {f)i, • • • , f)o — 1} inductively using lemma 5.3. 
By using the fact that the support of gh,uj is compact, we find that no longer 

depends on h as soon as /i ^ /ik — 2, i.e., ^ h^ — 2. Moreover, if 

/i ^ /ik — 2, the iterative equation for s^^^(k) (5.31) simplifies into 


^ai,a2 




a' 


(7.40) 


We can therefore write out (5.31) quite explicitly; for l)i ^ /i ^ /ik — 2 

sW(k) = 


+ 


1 - S'/Ik,^^2 




9hi,-l,UJ [t-Wr 


^(hk—1) ^ 


9h\^,u) 


( hk-2 

h'=h 

(jjh]^,ui + 9h]^—l,Lii ~ 9hi^ —1,01^2 9h^,S^ 

(7.41) 

where all the functions in the right side are evaluated at k. Note that in order to get 
the two-point function defined in section 1, we must integrate down to /i = f)/?: S 2 (k) = 
5(t/3)(k). This requires an analysis of the second and third regimes (see sections 8.3 
and 9.3 below). We thus find 

'S 2 (k) = (fffe^,^(k) |ft^-i,a;(k)) (1 - cr(k) - o-</j^(k)) (7.42) 



where 


a(k) ;= Wf- Vkf ^^9h^k 


(7.43) 


and 


+ 9h^-i,Lo) ^ I 

/ftk-2 \ 

E 'b'"'’ 


II 


“ 1 “ “ T/TT^^^k — “ 

9h-^,u) + 9h-^—l,u) ~ 9h]^ — l,u)^2,(jj 9h\^,Li 

(7.44) 

/V (U'\ — — 

in which IT 2 with h' G { 1)2 + 1, • • • , 1)2 — 1} U {l)i -|-1, • • • , l)i — 1} should be interpreted 
as 0. 


2 - Bounding the error terms. We then use (7.1), (7.10) as well as the bound 
|(5/ik,a; +5/ik-l,a;)“^| ^ (const.) 2'^'^ (7.45) 

which follows from (7.29) and (7.27), in order to bound cj(k): 

|(T(k)| ^ (const.) 2^'^\U\. (7.46) 


73 



Furthermore, if we assume that 


E (if''(k) 


h'=i,p 


^ (const.) 2^^"|C/| 


(7.47) 


which will be proved when studying the second and third regimes (8.42) and (9.63), then 


cj<ft^(k)| ^ (const.) 2'"'^\U\. 


(7.48) 


3 - Dominant part of the dressed propagators. Furthermore, it follows 
from (7.32) that 


5/ik,a;(k) +5/ik-l,a;(k) = “ 


koko + Id" 


/ -ih o_ q r \ 

0 —iko ^ 0 

0 C -iko 0 

V C 0 0 -iko I 


(1 + a') (7.49) 


where we recall (7.30) 


ko ■■= Zh^ko, ko := Zh^ko, i := Vh^i 


(7.50) 


in which Zh^, zy^^ and Vh^ were defined in (7.26) and satisfy (see (7.27)) 

where c[^\ and are constants (independent of hk, U and e). Finally the error 
term a' is bounded using (7.38) and (7.34) 

|a'(k)| ^ (const.) ((1 + \U\\K\)2^'^ + 2^^-’^^). (7.51) 


4 - Proof of Theorem 1.1. We now conclude the proof of Theorem 1.1, under 
the assumption (7.47): we define 

zi := zi,^, zi:=ztj^, vi := v,,^ 

and use (7.24) to bound 

k/ik “^il ^ (const.) \U\2^'^, \zh^^ -zil ^ (const.) \U\2'^'^, \vh^, -t^il ^ (const.) \U\2^^ 

which we inject into (7.49), which, in turn, combined with (7.42), (7.46), (7.48) and (7.51) 
yields (1.14). 
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7.4 Intermediate regime: first to second 

1 - Integration over the intermediate regime. The integration over the 
intermediate regime between scales f)i and f)i can be performed in a way that is entirely 
analogous to that in the bulk of the first regime, with the difference that it is performed 
in a single step. The outcome is that, in particular, the effective potential on scale f}i 
satisfies an estimate analogous to (7.10) (details are left to the reader): 

y rfx x’“W'2*l{x) < 

. y dx |(x - < Ci2-5-”'|f/| (7,52) 

^ y* |(x-X2,)"*»'4>,„(x)| < 
for / ^ 3 and m ^ 3. 

2 - Improved estimate on inter-layer terms. In order to treat the second 
regime, we will need an improved estimate on 

J djC (7-53) 

where {i,j) are in different layers, i.e. {a, a') G {a,b} x {d,b} or {a, a') G {d,b} x 

{a,b}, h' ^ 1 ) 1 . Note that since is proportional to any 

contribution to ^'(o a') contain at least one propagator between different layers, 

i.e. g(h",u),(a,a') with h' < h" ^ f}o or g(h"),{a,a') with h" ^ 0, and (d, a') G {a, b} x 

{d,b} U {a, 6} X {a, 6}. This can be easily proved using the fact that if the inter-layer 

hoppings were neglected (i.e. 7 i = 73 = 0 ), then the system would be symmetric under 

'fk,a ^ V’k.a, V’k,b V’k^ '4’k,b ^ V’k.fe 

(h') 

which would imply that Wf , ,, = 0. The presence of at least one propagator between 

different layers allows us to obtain a dimensional gain, induced by an improved estimate 
on each such propagator. To prove an improved estimate on the inter-layer propagator, 
let us start by considering the bare one, g^h",u)),{a,a') with (d, a') G {a, b} x {a, b}U{d, b} x 
{a, b} and h' < h” ^ f)o (similar considerations are valid for the ultraviolet counterpart); 
using the explicit expression (2.17) it is straightforward to check that it is bounded as in 
(4.50), (4.49), times an extra factor e2~^ . We now proceed as in section 7.1 and prove 
by induction that the same dimensional gain is associated with the dressed propagator 
g(h",uj),(a,a')j with (d, d') G {a, 6 } x {d,b} U {d,b} x {a,b}, and, therefore, with (7.53) 
itself. 

2-1 - Trees with a single endpoint. We first consider the contributions 
^ 2 ^J{aa') trees T G with a single endpoint. The dhini) factor 
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in the estimate (5.23) can be removed for these contributions using the fact that they 
have an empty spanning tree (i.e. T(t) = 0), which implies that the z^^’s in the right 
side of (5.22) are all z^^^’s and not z^’s, and can be estimated dimensionally by a con¬ 
stant instead of dhiin)- Therefore, combining this fact with the gain associated to the 
propagator, we find that for all m ^ 3, 


/ 


dx 



{h') 

2,0;,(a, 



^ (const.) e2^ \U\. 


(7.54) 


2-2 - Trees with at least two endpoints. We now consider the contribu- 

(h') (h') 

tions IBo , to from trees with ^ 2 endpoints. Let vZ be the node that 

has at least two children that is closest to the root and let h* be its scale. Repeating 
the reasoning leading to (7.9), and using the fact that the x™ falls on a node on scale 
^ /i*, we find 


dx 


X 




0 

E : 

h*=h'-\-l 


t—mhl 




for any 6 £ ( 0 , 1 ), so that 


dx 




(h') 
2,u,{a,a') 


,jx) ^ (const.) e2®T'+mm(0,l-m)ft'|f;|2 (7 55 ) 


where 9' := 29 — 1 > 1. 


Combining (7.54) and (7.55), and repeating the argument in section 7.2, we con¬ 
clude the proof of the desired improvement on the estimate of g, and that 


dx 




(h') 


2,a;,(Q:,Q:') 


:x 


^ (const.) e 2 ^'^|C/|(l + 2 ““(°’i-”^)'*|C/|) 


(7.56) 


for m ^ 3. 


8 Second regime 

We now perform the multiscale integration in the second regime. As in the first 
regime, we shall inductively prove that gh,ui satisfies the same estimate as gh,u) (see (4.53) 
and (4.51)): for all m ^ 3, 

’ J dx \x^°x^^gh,U^)\ ^ (const.) 

^ l^/i,a;(k)| ^ (const.) 2’^+^^ 
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which in terms of the hypotheses of lemma 5.2 means 

Cfc = 2, cg = 1, dh{mo,mi,m2) = 2-^oh-{rm+m,)^ ^ 

Cg = (const.) and Cq = (const.) 2^'. 

Remark: As can be seen from (3.19), different components of gh,w scale in different 
ways. In order to highlight this fact, we call the {a, b} components massive and the 
{a,b} components massless. It follows from (3.19) that the Li norm of the massive- 
massive sub-block of gh,uj{^) is bounded by (const.) 2“^" (instead of 2“^, compare with 
(8.1)) and that the massive-massless sub-blocks are bounded by (const.) In 

the following, in order to simplify the discussion, we will ignore these improvements, 
even though the bounds we will thus derive for the non-local corrections may not be 
optimal. 

In addition, in order to apply lemma 5.2, we have to ensure that hypothesis (5.12) 
is satisfied, so we will also prove a bound on the 4-field kernels by induction (.^o = 3 in 
this regime, so (5.12) must be satisfied by the 4-field kernels): for all m ^ 3, 

|(x-X4)-kFf4(x)| ^ C'^lUldhim) (8.2) 

where C'^ is a constant that will be defined below. Note that in this regime, 

io = 3 > = 2 

Cfc Cg 

as desired. 


8.1 Power counting in the second regime 

1 - Power counting estimate. It follows from lemma 5.2 and (7.52) that for 
all m ^ 3 and some ci, C 2 > 0, 


1 


dx 




N=l^^n-(h) I 




'' -T 


pG^(t) pG^{r) 


\Pvq\—21 
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where li „>2 is equal to 1 if /^ > 2 and 0 otherwise, and O' ■.= 29 — \ > 0, so that 
1 


/3|A| 


I dx |(x-X 2 ,rB« „(x)| 


n 2(2-^) Yl 4-|C/|'--i2i 

bS 2 J(t) tiSl£(r) 


lv>2S he 




N 


Per. 


(h) 

T,l_r ,3 


I-Pun 1 — 2 / 


(8.3) 


2 - Bounding the sum on trees. By repeating the computation that leads 
to (6.11), noticing that if io = 3, then for v £ ill(T) we have 2 — |P„|/2 ^ —|P„|/ 6 , we 
bound 

E E n 2-^^cr (8.4) 

|PuQ 1=2/ 

for some constant C 3 > 0. Thus (8.3) becomes 
1 


/3|A| 


I dx |(x-X 2 ,rBW „(x)| < 2 '‘l"-'ISl(m) 2 -l'-‘''“ 

■E E 

EiLl(h-l)>/—l + 5jv,l 


,Nh, 


(C4|t/|)^*=i('‘“^^ 


(8.5) 


for some C 4 > 0. Note that, if ^ = 2, the contribution with iV = 1 to the left side admits 
an improved bound of the form 04^/1 (m) 2 ^ which is better than the corresponding 

term in the right side of (8.5). This implies 


and 


I ■ix|x’"iq 2 „(x)| < «2'-+'-5 i(m)|C/| 

' |(x-X4)’"B™„(x)|<C5&(m)(2'-+2'’''‘)|(7| 


(8.6) 


|(x- X2,)”'Bg^„(x)| < 2l'“+'“>‘"-'lffi(m)(c5|C/|)'-‘ 
for some C 5 > 0, with 1^3. By summing the previous two inequalities, we find 

^ y dx (x - X4)'"lTiE(^) ^ Cf,dh{m)\U\{l + C6\U\{e{he - h) + e®')) 


(8.7) 


I /3|A| 


dx 


(x - xnrw!d;ijx) <i 2''‘+'-)(2-o 


(8.8) 
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for some ce > 0 , which, in particular, recalling that in the second regime — h ^ 
—2he + C, for some constant C independent of e, implies ( 8 . 2 ) with 

C'f,-=Cf,{l + C 7 sup |C/|(e|loge|+e®')) 

\U\<Uo,€<eo 

for some C 7 > 0 . 

Remark: The estimates (8.3) and (8.4) imply the convergence of the tree expansion 
(5.8), thus providing a convergent expansion of ^ in U. 

Remark: The first of ( 8 . 8 ) exhibits a tendency to grow logarithmically in 2“^. This 
is not an artifact of the bounding procedure: indeed the second-order flow, computed 
in [ValO], exhibits the same logarithmic growth. However, the presence of the e factor 
in front of {he — h) ^ 2 | loge| ensures this growth is benign: it is cut off before it has a 
chance to be realized. 


8.2 The dressed propagator 

We now turn to the inductive proof of (8.1). We recall that (see (4.18)) 

gh,M = fh,u{^)A^l,{^) (8.9) 

where 

fji to 

4,(k) :=i(k) + /^;,,,(k)wf)(k)+ Wf')(k)+ Wf'^(k). 

h'=h-{-l 


We will separate the local part of A from the remainder by using the localization operator 
defined in (7.15) (see the remark at the end of this section for an explanation of why we 
can choose the same localization operator as in the hrst regime even though the scaling 
dimension is different) and rewrite (8.9) as 

gh,io{^) = (^Tift,a;(k)^ (1 + nAh^W^) (Lg[/j]^^(k))^ (8.10) 

where Tgj/j] is a shorthand for 

Lg[ft],a;(k) := {f^h+l,M - /^ft- 2 ,a;(k)) (£ife,^(k)) . 

Similarly to the first regime, we now compute CAh^u)A^) and bound 77H/i_^(k)Lg[/j] ,^(k). 
We hrst write the local part of the non-interacting contribution: 


£i(k) 


/ iko 

71 

0 

r \ 

71 

iko 


0 

0 

e 

iko 

73^ 

V ^ 

0 

73 r 

iko / 


( 8 . 11 ) 
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where 


( 8 . 12 ) 




1 - Local part. The symmetries discussed in the first regime (see (7.21) 
and (7.22)) still hold in this regime, so that (7.23) still holds: 


£lTf)(k' + p‘^,o) = - 


/ Kh'ko 0 Vh'C" \ 

lifih' Kh'ko Uh'i 0 

0 Vh'K Kh'ko 1‘iVh'i 

\ Vh'^ 0 iKh'K Kh'ko j 


(8.13) 


with Kh'ili'h'Kh'Kh'■,Vh'') £ Furthermore, it follows from (8.6) that if h' ^ f)i, then 

(8.14) 


iCh'l ^ (const.) \U\2^^, |4'| ^ (const.) \U\2^% \jlh'\ ^ (const.) |17|2'*', 

\vh'\ ^ (const.) \U\2*^^^, \Dhi\ ^ (const.) \U\2^~^. 


If t)i ^ h' ^ f}o, then it follows from (7.10) that 

\Ch'\ ^ (const.) \U\2’^', ICh'l ^ (const.) \U\2'^', \uh>\ ^ (const.) \U\2'^', 

and from (7.56) that 

Ifj'h'l ^ (const.) 2®^ |?7|, {Dh'l ^ (const.) 2^^ \U\ 
for some 6*' G (0,1). Injecting (8.11) and (8.13) into (4.18), we find that 


(8.15) 


(8.16) 


+ Pf,o) ~ ~ 


/ izhko -iirhh 0 vK* \ 

'yi'ihh izhko vK 0 

0 vK* izhko IsKC 

\ vK 0 737/1?* izhko J 


(8.17) 


where 


to to to 

4 := 1 + ^ Ch', dr/j := 1 + ^ fih', 7/^ ;= 1 + ^ K', 


h'=h 


h'=h 


h'=h 


to to 

Zh-='^+'^ Ch', t’/i := 1 + ^ ^h' 


(8.18) 


h'=h 


h'=h 


in which ?/j/, jlfi', 7^/, and Uh' with /i'G{f)i + l,-- - ,f)i — 1} are to be interpreted as 
0. By injecting (8.14) through (8.16) into (8.18), we find 


\mh - 1| ^ (const.) \U\, \zh - 1| ^ (const.) \U\, \zh - 1| ^ (const.) \U\, 
\vh - 1| ^ (const.) |f7|, \vh - 1| ^ (const.) \U\. 


(8.19) 
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2 - Dominant part of Furthermore, we notice that the terms propor¬ 

tional to Zh or Vh are sub-dominant: 


+ Pf,o) ~ ^^h,LoO^' + Pf,o)(^ + ^3(k0) 


( 8 . 20 ) 


where 


+ Pf,o) ~ ~ 


/ 

0 

lirhh 

0 

VhC \ 


Hrhh 

0 

Vh^ 

0 


0 

VhC 

izhko 

0 

V 

vhi 

0 

0 

izhko / 


( 8 . 21 ) 


Before bounding < 73 , we compute the inverse of (8.21), which is elementary once it is put 
in block-diagonal form: using proposition C.l (see appendix C), we find that if we define 


7i:=m/i7i, ko := Zhko, ^ := VhC 


then 




1 _ / 1 Mh(k)^ \ / a 




0 -at\k) 


where 


H ■=- 


a}^) - 

and 


0 7i ^ 

7f^ 0 


^(Pf,o + 1^0 •— 3 


71 


ifkl + 


1 0 
Mh{k) 1 


iiiko (rf 

ij^ko 




In particular, this implies that 

1-1 I'U' 


+p‘f,o)I ^ (const. 


2-h, 

h + h(^ 
2 —^ 


_ h+h^ 

2 ^ 

o—h 


( 8 . 22 ) 

(8.23) 

(8.24) 

(8.25) 

(8.26) 


in which the bound should be understood as follows: the upper-left element in (8.26) is 
the bound on the upper-left 2 x 2 block of and similarly for the upper-right, 

lower-left and lower-right. In turn, (8.26) implies 


, . , . . / h — he 3he — h\ 

|(T3(k )| ^ (const.) (^2 2 +2 2 j. 


(8.27) 


3 - Irrelevant part. The irrelevant part is bounded in the same way as in the 
first regime (see (7.36)): using (8.17) and the bounds (8.14) through (8.16), we find that 
for m ^ 3 and f )2 ^ ^ h' ^ f)i. 


dx 




^ 2 '"'^ 5^/1 (m)(const.) |17| 


(8.28) 
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(8.29) 


and for ^2 ^ h ^ ^ h' ^ {)o, 

J dx *hg[h],U^) ^ 2^-=5'fe(m)(const.) \U\. 

Therefore, using the fact that 



x™7^(5^ >i=Lg[/j]_^(x)| ^ 2^'= S'/I(m)(const.) 


we find 

j dx\x"^nAh,^ >i=Lg[/,]^^(x)| ^ 2^'=S/i(zZl)(const.) (1 + \h\\U\). 


(8.30) 


4 - Conclusion of the proof. The proof of (8.1) is then concluded by inject¬ 
ing (8.21), (8.27), (8.20) and (8.30) into (8.10). 

Remark: By following the rationale explained in the remark following (7.17), one may 
notice that the “correct” localization operator in the second regime is different from 
that in the first. Indeed, in the second regime, {k — ppQ)dk scales like 22 !^“^ ) instead of 
2 h-h first. This implies that the remainder of the first order Taylor expansion of 

1 T 2 ^ ^ is bounded by 2^ instead of in the first regime, and is therefore marginal. 

However, this is not a problem in this case since the effect of the “marginality” of the 
remainder is to produce the \h\ factor in (8.30), which, since the second regime is cut 
off at scale She and the integration over the super-renormalizable first regime produced 
an extra (see (8.30)), is of little consequence. If one were to do things “right”, one 
would define the localization operator for the massless fields as the Taylor expansion 
to second order in k and first order in ko, and find that the \h\ factor in (8.30) can be 
dropped. We have not taken this approach here, since it complicates the definition of 
C (which would differ between massive and massless blocks) as well as the symmetry 
discussion that we used in (8.17). 


8.3 Two-point Schwinger function 

We now compute the dominant part of the two-point Schwinger function for k 
well inside the second regime, i.e. 


fji-i 

k G := IJ suppfh,ej 

/l=f )2 + l 


Let 


hk := max{/r : fh,ui0A) / 0} 


so that if /i 0 {/ik, /ik — 1}) then fh,u)0^) = 0. 
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1 - Schwinger function in terms of dressed propagators. Recall that the 
two-point Schwinger function can be computed in terms of the effective source term 
defined in (5.27), see the comment after lemma 5.3. Since /ik ^ f)i, X^^^ is left invariant 
by the integration over the ultraviolet and the first regime, in the sense that X^^^^ = 
withAf('jo) given by (7.39). We can therefore compute X^^'^ for h £ {f) 2 , • • • , f)i —1} 
inductively using lemma 5.3, and find, similarly to (7.42), that 

S 2 (k) = + |ftk-l,a;(k)) (1 - cr(k) - fT<hk(k)) (8.31) 

where 

u(k) := IRf + {gh^,u + 

(8.32) 

and 


\h'=hf3 


~ I ~ ~ T/TT'^^k —~ 

I — — ^2 cj 9h]^,o 


(8.33) 


2 - Bounding the error terms. We now bound cj(k) and (T<ft^(k). We first 
note that 

^ (const.) (8.34) 

which can be proved as follows: using (8.9), we write gh^,u) = Ak^hkoj 

lftk-l,aJ = Ak-l^fck,1.^(1 + 

Therefore, noting that //i,^(k) -|- fh^_^(k) = 1, we obtain 
{ghi^,u}+ghi^-i,Lo) = Ak 1 +Ak-i((k + /^/ik-i^2^*' 

Now, by (8.6), we see that |lT2^''~^^(k)^/0^^t^(k)| ^ (const.) 2^", which implies (8.34). 
By inserting (8.34), (8.6) and (8.1) into (8.32), we obtain 

|cr(k)| ^ (const.) 2^^\U\. (8.36) 

Moreover, if we assume that 

h 

^ lT 2 ^'''^(k) ^ (const.) 2^^'17/1 (8.37) 

^' = f)/3 
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4 - Proof of Theorem 1.2. We now conclude the proof of Theorem 1.2, under 
the assumption (8.37). We define 

:= (1 +Cr'(k)) +|hk-l,a;(k))"^ 

(i.e. the inverse of the matrix on the right side of (8.39), whose explicit expression is 
similar to the right side of (8.21)), and 

m2:=m(,2, Z2 := Zf,^, V2 := 

and use (8.14) to bound 

\mh^^ - m 2 \ ^ (const.) \U\2^'^, \zh^^ - Z 2 I ^ (const.) |C/||/i,|2''% 

-^21 ^ (const.) |C/|25('*k+/id 

so that 

{B^,{k) - ^ (const.) \U\\h\2’^^ 

which implies 

B-^{k) = B-^\k){l + 0{\U\\h,\2^^)). (8.41) 

We inject (8.41) into (8.39), which we then combine with (8.31), (8.36), (8.38) and (8.40), 
and find an expression for S 2 which is similar to the right side of (8.39) but with 
replaced by f) 2 - This concludes the proof of (1.18). Furthermore, the estimate (1.23) 
follows from (8.19), which concludes the proof of Theorem 1.2. 
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5 - Partial proof of (7.47). Before moving on to the third regime, we bound 
part of the sum on the left side of (7.47), which we recall was assumed to be true to 
prove (1.14) (see section 7.3). It follows from ( 8 . 6 ) that 




h'=i)2 


^ (const.) 2^^-\U\. 


(8.42) 


8.4 Intermediate regime: second to third 

In the intermediate regime, we integrate over the first scales for which the effect 
of the extra Fermi points cannot be neglected. As a consequence, the local part 

of A(, 2 ,aj(k) is not dominant, so that the proof of the inductive assumption (8.1) for 
/i = f )2 must be discussed anew. In addition, we will see that dressing the propagator 
throughout the integrations over the first and second regimes will have shifted the Fermi 
points away from p*^^. by a small amount. Such an effect has not been seen so far because 
the position of p'^Q is fixed by symmetry. 

1 - Power counting estimate. We first prove that 

j dx |x”' 5 (, 2 ,^(x)| ^ (const.) 2 “^ 2 ^(, 2 (m). (8.43) 

The proof is slightly different from the proof in section 8.2: instead of splitting 
according to ( 8 . 10 ), we rewrite it as 

(i(k) + £*m„ 2 ,-(k))~' (l + (L0[^2],-(k)))~' (8.44) 

(this decomposition suggests that the dominant part of A[, 2 ,a; is A + CWt, 2 ,u) instead of 
>CA(, 2 ,aj) in which we recall that A = 

im^ 2 ,a.(k) := if, 2 ,a;(k) - i(k) 


and 

L 0 [f, 2 ],<..(k) := j /^f, 2 +l,a;(k) - ^ /^(, 2 _ 2 ,a;,i(k) j (i(k) +£fmt, 2 ,^(k)) (k). 

\ ie{0,i,2,3} / 

We want to estimate the behavior of (8.44) in which we recall is a ball with 

four holes around each p%j, j = 0,1,2,3. The splitting in (8.44) is convenient in 

that it is easy to see that A(k) + £ 2 Ilt, 2 ,a;(k) satisfies the same estimates as A(k); in 
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particular, via proposition B.l (see appendix B), we see that det (^(k) + £ 2 ir(, 2 ,cj(k)) ^ 
det^(k) • (1 + 0{U)) on so that for all n ^ 7 and k G 

(i(k) + £im„ 2 ,-(k))~' ^ (const.) (8-45) 

and, moreover, for m ^ 3, 

j dx |x™772B(,2,^ *L 0 [^ 2 ],i^W| ^ (const.) \U\\he\2^‘^i,^{rn). (8.46) 

The proof of (8.43) is then concluded by injecting (8.45) and (8.46) into (8.44). We can 
then use the discussion in section 8.1 to bound 

I dx \x^W^^l{x)\^C22^-+>^^d-^^{m)\U\ 

< ^ y dx [x-Xi)'^wQjx) i^C 2 ^^{m)\U\ (8.47) 

dx |(x-X2z)”^1t£Jx)| ^ 2«^+'^^)(2-05^^(^)(C2|[/|)'-i 

for some constant (72 > 1 . 


2 - Shift in the Fermi points. We now discuss the shift of the Fermi points, 
and show that 5 ^( 12 ,cj has at least 8 singularities: p'^g and for j G {1,2,3} where 


_l_ (o,0,a;A[,2) 


(8.48) 


_ rp-u)~{0J,t)2) ~i‘^A2) _ ™ - (<^,l)2) 

>^F,2 ~ ^ nF,! ’ nF,3 ~ t^F,! 

in which denotes the spatial rotation by ±2 f/ 3; and that 


(8.49) 


|A(, 2 | ^ (const.) e‘^\U\ 


(8.50) 


(note that (8.49) follows immediately from the rotation symmetry (2.33), so we can 
restrict our discussion to j = 1 ). 


Remark: Actually, we could prove in this section that g^t) 2 ,Lo has exactly 8 singularities, 
but this fact follows automatically from the discussion in section 9, for the same reason 
that the proof that the splittings (7.18) and (8.10) are well defined in the first and second 
regimes implies that no additional singularity can appear in those regimes. Since the 
third regime extends to /i —)■ — 00 , proving that the splitting (8.10) is well defined in the 
third regime will imply that there are 8 Fermi points. 


86 



We will be looking for in the form (8.48). In particular, its ko component 

vanishes, so that, by corollary B.2 (see appendix B), Afjj solves 








In order to solve (8.51), we can use a Newton iteration, so we expand around 0; 

it follows from the symmetries (2.35) and (2.36) that 

D^2,Ul{^t)2) ~ ^i)2 +^^t2^t2 (8.52) 

with (M[, 2 ,yt) 2 ) G independent of cj. Furthermore by injecting (7.10) and (8.6) 
into (8.51), we find that 


^2 = ^7173 + 0{^\U\) + O(e^), Mf ,2 = 0{e^\U\) (8.53) 


and 




^ (const.) . 


(8.54) 


Therefore, by using a Newton scheme, one finds a root Af ,2 of (8.51) and, by (8.53) 
and (8.54), 

|A[, 2 | ^ (const.) e^\U\. (8.55) 

This concludes the proof of (8.48) and (8.50). 


9 Third regime 

Finally, we perform the multiscale integration in the third regime. Similarly to the 
first and second regimes, we prove by induction that gh,ui,j satisfies the same estimate 
as gh,uj,j (see (4.56) and (4.54)); for all m ^ 3, 

’ J dx ^ (const.) 2->^-^oh-7n,^{h-h,) 

i ll/i,a;,i(k)| ^ (const.) ^ ^ 

which in terms of the hypotheses of lemma 5.2 means 

Ck = 3, c, = 1, dh{mo,mi,m2) = 2-^oh-i^r+^2)ih-h.)^ 

Cg = (const.) and Cq = (const.) 2“^^'. 
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Remark: As in the second regime, the estimates (9.1) are not optimal because the 
massive components scale differently from the massless ones. 


Like in the first regime, 


4 = 2 > 



3 

2' 


9.1 Power counting in the third regime 

1 - Power counting estimate. By lemma 5.2 and (8.47), we find that for all 


m ^ 3 
1 


/3|A| 


dx 


E E E 

It- pp-p(^) 

\Pvq\=21 

(^l22h.)iV-l 

n 2(^ |A;|) (C22 I 1)2(^L l)^e 


SO that 

1 I 


/3|A| 


dx 


1.x - X2« J 




y~^ yy yy JJ gL|[/|max(i,z„-i)^ 

PG^(r) i;G(£(r) 

|Pt,Q 1=2/ 

(9.2) 


A'" —1 g'T-('») L pc'p(^) 


2 - Bounding the sum of trees. We then bound the sum over trees as in 
the first regime (see (7.4) and (7.8)): if f ^ 2 then for 9 G (0,1) and recalling that 
f )2 = 3/ig + const. 


N 

^ JJCpL (9.3) 

\Pvq\=21 

and if / = 1 then 

N 

Y Y n (9-4) 

I I =2 
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Therefore, proceeding as in the proof of (7.10) and (7.12) we find that 
[ dx x™lTi2_._„(x) ^ 2^(^-^^^:Sh{rn)Ci\U\ 


(9.5) 


and 


f (x-x 4 )™it'S.-„(x) ^ dh{m)Ci\u\ 


/3|A| J 1^— "4,^,j,o 

rW 


[ ^/dx (x - X 2 z)™ 1 T 4 ';^„(x) ^2(3-2O/^+20(/^-3/^d + ( 2 ^-l)/^. 5 ^(^)(C'l|[/|)^-l 


(9.6) 


for / > 3 and m ^ 3. 


Remark: The estimates (9.2), (9.3) and (9.4) imply the convergence of the tree expan¬ 
sion (5.8), thus providing a convergent expansion of ^ in U. 


9.2 The dressed propagator 

We now prove (9.1). We recall that (see (4.23)) 

gh,u,j (k) = fh,L,,j (k)ift,L,i (k) (9-7) 

where 

t2 

4^j(k) = i(k) + ^ wPik) 

h'=h+l 

ti 1)0 

+ 5: af\u) + 5: 

h'=t)2 /i'=f)i 

1 - j = 0 case. We first study the j = 0 case, which is similar to the discussion 
in the second regime. We use the localization operator defined in (7.15) and split gh,uj,o 
in the same way as in (8.10). We then compute CW 2 and bound TZA^^uj 


1-1 - Local part. The symmetry considerations of the first and second 
regime still hold (see (7.21) and (7.22)) so that (7.23) still holds: 


£wi'^')(k' + p‘^,o) 


/ Kh'ko -fifth' 0 Uh'f* \ 

lifh' Kh'ko Uh'f 0 

0 Vh'f* Kh'ko 1‘iVh'i 

\ Vh'f 0 -fKh'K Kh'ko / 


(9.8) 


with Kh'i fh'-, kh', Ch', ^h’) £ The estimates (8.14) through (8.16) hold, and it follows 
from (9.5) that if /i' ^ 1)2, then 


\Ch'\ ^ (const.) |[/|2'^'-2^% \Ch'\ ^ (const.) |[/|2'^'-2^% \fih'\ ^ (const.) 

^ (const.) \U\2^'-^% ^ (const.) |C/|2^'-2'^^ 

(9.9) 
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Therefore 


+ P^^o) ~ ~ 


/ izhko 7 im/i 0 VhC \ 

-fifhh izhko VhC 0 

0 VhC izhko -f3Vh^ 

\ Vhi 0 ^3VhC izhko J 


(9.10) 


where Zh, Zh, nih, Vh and Vh are defined as in (8.18). and are bounded as in (8.19): 

\mh — 1| ^ (const.) \U\, \zfi — 1| ^ (const.) 1^71, \zh — 1| ^ (const.) \U\, 

\vh - 1| ^ (const.) \U\, \vh - 1| ^ (const.) \U\. 


(9.11) 


1-2 - Dominant part of Furthermore, we notice that the terms 

proportional to Zh are sub-dominant; 


71Ah,aj,o(k -|- Pf^o) ~ -|- p‘pq)( 1 -|- CT4(k )) 


(9.12) 


'2^/i,aj,o(k^ + Pf,o) ~ 


(9.13) 


where 

/ 0 'yifhh 0 VhC \ 

'yirhh 0 Vh^ 0 

0 VhC izhko 

\ Vhi 0 -fzVhC izhko J 

Before bounding 1 T 4 , we compute the inverse of (9.13) by block-diagonalizing it using 
proposition C.l (see appendix C): if we define 


ko := Zhko, 71 := rfih'yi, ^ := Vh^, C ■= 


(9.14) 


then for k G B 
^'C^h,i.j,o(k)^ 


(h,LLj,0) 

kL ’ 
-1 / 


1 

0 1 


“/i ,0 ^ 

-(m) 

%,0 


0 47 (k); V ^M(k) 1 


1 0 


(1 + 0 ( 2 “')) 

(9.15) 


where 


- - 

“/i,o 


-1 


0 7i 

7r' 0 


’ “m (P+o + k') -77 


-iko 73^ 


4 ' + 7 |llPV 73 e -iko 


(9.16) 


(the 0 ( 2 ^ ^^') term comes from the terms we neglected from that are of order 

2-3+-) 


In 


Mh,o{pk,o + k') Q ^ 

particular, this implies that, if (k' + p'^q) G then 

- / ‘2~2he 

l'^^hi,o(k' + PF,o)l ^ (const.) 2 + 


(9.17) 


(9.18) 
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in which the bound should be understood as follows; the upper-left element in (9.18) is 
the bound on the upper-left 2x2 block of and similarly for the upper-right, 

lower-left and lower-right. In turn, (9.18) implies 

|cj 4 (k')| ^ (const.) (9.19) 


(h') - 

1-3 - Irrelevant part. We now bound 'R-W^Jq * o] ™ the same way 


as in the second regime, and find that for m ^ 3, if h ^ /i' ^ f)o, then 






^ 2 


h—2hf 


5^/j(m)(const.) \U\ 


so that 

j dx\x^nAh,cj,o*^Q[h],Lo,oi^)\ him) (const.) {l + \h\\U\). 

This concludes the proof of (9.1) for j = 0. 


(9.20) 


(9.21) 


2 - J = 1 case. We now turn to the case j = 1 (j = 2,3 will then follow by 
using the 27r/3-rotation symmetry). Again, we split gh,u),i in the same way as in (8.10), 
then we compute CW 2 ^ and bound 'RAh^u,i^Q[h\,ui,i- Before computing CAh^u),i and 
bounding 'R-Ah^u),i^Q[h],u),i^ we first discuss the shift in the Fermi points (i.e., the 

singularities of ^(k) in the vicinity of p^j^^), due to the renormalization group flow. 


2-1 - Shift in the Fermi points. We compute the position of the shifted 
Fermi points in the form 

p^-f) = p-i + (0,0,a;A;,) (9.22) 

and show that 

|A/i| ^ (const.) e^\U\. (9.23) 

The proof goes along the same lines as that in section 8.4. 


Similarly to (8.51), A/^ is a solution of 


Dh,ui,liA>.h) •— ^ft,aj,l,(fe,a)(PF,i (Pf,! ’)Ah,u),l,(b,d)iPF^[ 




U 


= 0 . 


(9.24) 


We expand Dh,uj,i around A/j+i: it follows from the symmetries (2.35) and (2.36) that 
Dh,uj,ii^h) = Mh + ojYh(Ah - Ah+i) + (Ah - Ah+i)'^R]^j^^^(Ah) (9.25) 

with (Mh,Yh) G independent of uo. Furthermore, 


Mh = Ilh,tj,i(A/i+i) = Dh,u},i(Ah+i) — Dh+i,uj,i(Ah+i) 
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so that, by injecting (9.5), (7.10) and (8.6) into (9.24) and using the symmetry structure 
of ^h^ij^i(k) (which imply, in particular, that |A/i^^^i(k)| ^ (const.) e in we 

find 

\Mh\ ^ (const.) (9.26) 


and 


as well as 


27173 + 0(e'|C/|) + 0(6"). 


^ (const.) {l + e\U\\h\). 


(9.27) 


(9.28) 


Therefore, by using a Newton scheme, we compute Ah satisfying (9.24) and, by (9.26), 
(9.27) and (9.28), 

\Ah - A/i+il ^ (const.) (9.29) 


This concludes the proof of (9.22) and (9.23). 


2-2 - Local part. We now compute CAh^uj^i- The computation is similar to 
the j = 0 case, though it is complicated slightly by the presence of constant terms in 
Ah,ui,i- Recall the x-space representation of Ah^u),i (4.42). The localization operator has 
the same definition as (7.15), but because of the shift by in the Fourier transform, 

its action in k-space becomes 

+ (k - p^f^) • 9kih,a,,l(P^f^)- 

In order to avoid confusion, we will denote the localization operator in k space around 
by Ch- 

2-2-1 - Non-interacting local part. As a preliminary step, we discuss 
the action of t on the undressed inverse propagator A(k). Let us first split A(k) into 
2x2 blocks: 

^ i00(k) ) 
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in terms of which 


+ p'-pf) = Ua*HK + pfe"’) 


iko 7i \ 
7i iko ) 


m 


( 0 ) 


m 


( 0 ) 


+ 


^i, 


x + <^wfh[. 


+ 


kl 

0 


+ ujw^^h[ 


ChA<^^ik[ + 


ik'^,1,0 + iiv^h^kl^ + ujw'j^’k[^y)) \ 

(9.30) 


,(0)i 


where 


= 7173 + C>(Afe), ^ 


4°^ = ^ + Afe), = ^ + 0(e^ A,,). 


(9.31) 


" " ^ (h') 

2-2-2 - Local part of 1V2- We now turn our attention to ChW^ ■ In 

^ ih') 

order to reduce the size of the coming equations, we split II 2 into 2x2 blocks: 


Wf=: 


W. 


(h'm 


w. 


(h')(<fi 


2 ^*'^2 


The symmetry structure around p^(^^ is slightly different from that around p'^g- 
deed (7.21) still holds, but p^(^^ is not invariant under rotations, so that (7.22) becomes 

Pp’i — Pf,i ~ ^vPpi — ^hPp^i — ^Ppi — • f9.32j 

It then follows from proposition F.l (see appendix F) that for all (/,/O £ 

^_( + ii^U'ik[x + ujwlf^kly) \ 

V + i-i^h(A,X + ^^i(lk'l,y) J 

(9.33) 
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(9.34) 


f ff t S f A 

with il-i-fj i,Ch'^h'I’^h' i) £ III addition, by using the parity symmetry, it follows 
from (F.IO) (see appendix F) that the block is equal to the block. Furthermore, 
it follows from (9.5) that for h' ^ f) 2 , 

^ (const.) |^//'| ^ (const.) 

Wh(i\ ^ (const.) \U\2^'-^% Iro^/'I ^ (const.) 

If i )2 ^ h' ^ 1)1, then it follows from (8.6) that 

\Ci(i\ ^ (const.) |C/|2^s 

\ulf[\ ^ (const.) ^ (const.) |I/|2i(^'+^^) 


and because 11 ^ 2 ^ ^(Pfo) ~ *-* ^cid \p'P(‘'' — P^ol ^ (const.) 2^'*% by expanding 1 ^ 2 '^ '* to 
first order around p'pg, we find that it follows from (8.6) that 

^ (const.) 25(^'+^^)22'*^|C/|. 

Finally, if fii ^ /i' ^ f)o, then it follows from (7.10) that 

\Ci(i\ ^ (const.) \U\2^', 

Wlfi\ ^ (const.) \U\2^', ^ (const.) |17|2^' 

^ (h') 

and by expanding IF 2 ^ to first order around p'^g, we find that 

U(,\ ^ (const.) |t/|2'^'+2'^^ 

By using the improved estimate (7.56), we can refine these estimates for the inter-layer 
components, thus finding: 

\lJ-{Jl\ ^ (const.) \U\2^h'+3hc^ 

^ (const.) |t/|2®^'+^% ^ (const.) |[/|2®^'+^% (9.39) 

ic^pl = \citi\ ^ (co^st.) \U\ 2 ^’^'+^^ 

for all / G 


.{uj,h) 


(9.35) 


.{h') 


(9.36) 


(9.37) 


(9.38) 


2-2-3 - Interacting local part. 

with (9.30), we find 


Therefore, putting (9.33) together 


(k'l+pfe"^) 



/ 

izfiko 

\ ^ti + Ki^i 

7lKi+<f) 

m,kt, 

izl*,ko 

izl^iko 

iz'l^iko 


izl*,ka 

73 (mS + Ktt 


izt)h 


(9.40) 
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with 




-'h,! 


h,l '^l,y 


for (/,/') e {(/),^}2, and 


(0) 

K:i — K 


llo 


1)0 


+ ™a-i + ^E4h 


.__( 0 ) 


f)o 


73 


h’=h 


71 


+ E '■i? 


1’ 




h'=h 


■'7=i + Eci/.. 


f)0 


i-p ._ 


f)0 


V 


h'=h 


h'=h 


1 l)o 

^h,l •- + ^3 ^'7 


h'=h 

to 


— - — 
71 


to 


E7i.. 


._ „(") _ 


to 

E 


/i',i’ 


1 vu 

<p<p ._ .,,(0) + Y" 
73 ^ 




^4>4> 


— 


h'=h 

^0 


wEi •= ~ ^h',1 


■'ft',!’ '^h,l 


^<f> (0) 

wi, := w 


’ '^hA 


h'=h 
to 

h + £ ^itr 


h'=h ' h’=h 

Furthermore, using the bounds (9.34) through (9.39), 

I’^h'i “ + l^-i^i “ 1| + I’^i'^i “ ^ (const.) e^\U\, 

\^h!i - M ^ (const.) \U\, ^ (const.) |loge|e|[/|, 

I^M " ^ (const.) \U\, 

+ l^i^il + l^h '^1 “ ^ (const.) \U\. 


h'=h 


(9.41) 


(9.42) 


2-2-4 - Dominant part of ChAh,oj,i‘ Finally, we notice that the terms 
in (9.40) that are proportional to or Kf^\ are subdominant: 

ChAh,+ P^F,l'^) = thAh,oj,liK + P^i^^)(l -F fT4,i(k'J) (9.43) 

where 

^hAh,u:,i{K +P^f^) 


/ 

0 

71 

0 

i4> I 

'''-h,! ^h,l 

\ 


7ini£ 

0 

™i* + Kf* 

0 



0 

?</< I 

'''-h,! ^h,l 

‘4*1=0 

73«. + <t) 


V 

_|_ 

0 


‘4*1=0 

/ 


(9.44) 
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Before bounding < 74 ^ 1 , we compute the inverse of (9.44) by block-diagonalizing it using 
proposition C.l (see appendix C): if we define 




ko := zf'jlko, Hi := xi := (9-45) 


then for k G . 




where 






_ 1 Mt (k) C 


0 1 


0 4™)(k) i V MhA\^) 1 


7173 


1 0 


(1 + 0(2^-^’^^)) 
(9.46) 


- - 


---1 ^0 )’ 4 ?(P ‘^4 + k'i) := 


_ iko 73 X 1 

7i " 9 J' -u ^2 _^^2|^^|2 73^^ ikQ 

(the term comes from the terms in of order and 

0 

71 V 9 H; 

In particular, this implies that if (k^ -|- p^4) G then 


+ ki) :— — — 1 p 

7 i \ U r .1 


I [hhAh,+ P^ 4 )] ^1 ^ (const.) 


<2^he — h 2^he — h 
nhe—h o—h 


(9.47) 


(9.48) 


(9.49) 


in which the bound should be understood as follows; the upper-left element in (9.49) is 
the bound on the upper-left 2 x 2 block of ^hAj^lj i; and similarly for the upper-right, 
lower-left and lower-right. In turn, using (9.49) we obtain 

|o' 4 ,i(k'i)| ^ (const.) e(l -f | loge||C/|). (9.50) 


2-3 - Irrelevant part. Finally, we are left with bounding TZAh,u,i^Q[h],uj,i^ 
which we show is small. The bound is identical to (9.21); indeed, it follows from (9.46) 
and (9.49) that for all m ^ 3, 

j dx |x™L 0 [;,]^^^i(x)| ^ (const.) 2-^^h{m) 

so that 

j dx\x^nAh,uj,i *Ih 0 [/i],,^,i(x)| ^ him) (const.) (I -f \h\\U\). (9.51) 
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3 - j = 2,3 cases. The cases with j = 2,3 follow from the 27r/3-rotation 
symmetry (2.33); 




1 


0 


0 Trp,, ~(lAJ,h) 

J- ^j+PF,j-uj 


[ 0 7-t 


T\c' 


where T and 7k were defined above (2.33), and 


(9.52) 


9.3 Two-point Schwinger function 

We now compute the dominant part of the two-point Schwinger function for k 
well inside the third regime, i.e. 

fj2-i 

y suppfh,u;,j- 

^=ll/3 + l 


Let 


hk := max{/r : fh,u,j{^) / 0} 


so that if /i 0 {/ik, /ik “ 1}) then fh,uj,j0^) = 0. 


1 - Schwinger function in terms of dressed propagators. Recall that 
the two-point Schwinger function can be computed in terms of the effective source term 
see (5.27) and comment after Lemma 5.3. Since /ik ^ ^ 2 , is left invariant 
by the integration over the ultraviolet, the first and the second regimes, in the sense 
that with given by (7.39). Therefore, we can compute for 

/i G {f)/ 3 , • • • ,1)2 — 1} inductively using lemma 5.3, and find, similarly to (7.42) and (8.31), 
that 

'S2(k) = {§hi,,Lo,j(k) + §hi,-i,co,j(k)) (1 - cr(k) - fT<hk(k)) (9.53) 

where 


(9.54) 

and 


h\,—2 


{h') 


C<hk(k) (^1 {ghy„U,j+ghy,-l,U,j) ^S'/IkI ^2 

\^'=ll/3 

- I - - Ti;-(^k~l)- \ 

9h-^,u),j + gh^-l,u),j — ) * 


Similarly to (8.34), we have 

|(7/ik,<.^j ^ (const.) 


(9.55) 

(9.56) 
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and, by (9.5) and (9.1), we have 


cj(k)| ^ (const.) 2^'^ 
a<,,Jk)| ^ (const.) 


(9.57) 


2 - Dominant part of the dressed propagators. We now compute gh^,uj,j + 


2-1 - j = 0 case. We first treat the case j = 0. It follows from (the analogue 
of) (8.10), (9.12) and (9.15), that 


\ 0 


*,Lo(k) 



1 

Mh„o(k) 


(1 + cro(k)) 


(9.58) 

where were defined in (9.17) and (9.16), and the error term cjq 

can be bounded using (9.21) and (9.19); 


ao(k)| ^ (const.) {2-^^ + \hk\\U\). (9.59) 


2-2 - j = 1 case. We now consider j = 1. It follows from (9.43) and (9.46) 


that 


5/iki‘*.’,o(k) + ff/ik—i,tiJ,o(k) — 




1 <^,i(k) Cd 


0 aS(k) ) V 1 


(l + aUk)) 

(9.60) 


where and were defined in (9.48) and (9.47), and the error term a[ 

can be bounded using (9.51) and (9.50): 

|a((k)| ^ (const.) ( 2 ^^(l + |/r,||t/|) + 2 '*k- 2 /^.( 2 -/*.+ |/jj^||[/|)^ . (g.ei) 

2-3 - j = 2,3 cases. The cases with j = 2,3 follow from the 27r/3-rotation 
symmetry (2.33) (see (9.52)). 


3 - Proof of Theorem 1.3. We now conclude the proof of Theorem 1.3. 
We focus our attention on j = 0,1 since the cases with j = 2,3 follow by symmetry. 
Similarly to section 8.3, we define 

^/Jkj(k) := (1 + f^j(k)) (0h^,a;j(k) + 
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(i.e. the inverse of the matrix on the right side of (9.58) for j = 0, (9.60) for j = 1, 
whose explicit expression is similar to the right side of (9.13) and (9.44)), and 

hi3,o := ^3,0 := ^^3,o := '^1)2) '^3,0 := h[,2, 

a - 4>4' i4> 4’4’ 

m 3 ,i := p 1113,1 := m 3 ,i := ^3,1 := 4 ;,i, 

— fd) - E(t> ~ ~ 

V3,1 ■■= W3,1 := r;3,i := %,i, w^3,i := 

and use (9.9) and (9.34) to bound 

\fhh^ - fh 3 fl\ + - ?fi 3 ,i| + - m 3 p| ^ (const.) |[/|22^k-3/i.^ 

1 “ "^3,il ^ (const.) |[/| 22 ^k- 2 /ie^ 

k/ik - ^3,i| + khk.i “ ^3,o| ^ (const.) 

\vh^ - U3,o| + - t’3,i| + - W3,i\ ^ (const.) \U\2^'-~^% 

|h/ik “ ^3,o| + - h3,i| + ~ ^3,il ^ (const.) |[/|2'"k-2/i. 

so that 

(B^ 2 o(k) -i?hk,,(k))il^-;,(k)| ^ (const.) \U\2’^--^^^ 

which implies 

+ O(|t/|2">^-2^0)- (9-62) 

We inject (9.62) into (9.58) and (9.60), which we then combine with (9.53), (9.57), (9.59) 
and (9.61), and hnd an expression for S 2 which is similar to the right side of (9.58) 
and (9.60) but with /ik replaced by 1)2- This concludes the proof of (1.24). Furthermore, 

the estimate (1.29) follows from (9.11) and (9.42) as well as (9.31) and (9.23), which 

concludes the proof of Theorem 1.3. 

4 - Proof of (7.47) and (8.37). In order to conclude the proofs of Theorems 1.1 
and 1.2 as well as the Main Theorem, we still have to bound the sums on the left side 
of (7.47) and of (8.37), which we recall were assumed to be true to prove (1.14) and (1.18) 
(see sections 7.3 and 8.3). It follows from (9.5) that 

fj2 

Wi^'\k) ^ (const.) 2^'^'|17|. (9.63) 

h'=t)p 

This, along with (8.42) concludes the proofs of (7.47) and (8.37), and thus concludes the 
proof of Theorems 1.1, 1.2 and 1.3 as well as the Main Theorem. 


10 Conclusion 

We considered a tight-binding model of bilayer graphene describing spin-less fermions 
hopping on two hexagonal layers in Bernal stacking, in the presence of short range inter- 


99 



actions. We assumed that only three hopping parameters are different from zero (those 
usually called 70,71 and 73 in the literature), in which case the Fermi surface at half¬ 
filling degenerates to a collection of 8 Fermi points. Under a smallness assumption on 
the interaction strength U and on the transverse hopping e, we proved by rigorous RG 
methods that the specific ground state energy and correlation functions in the thermo¬ 
dynamic limit are analytic in U, uniformly in e. Our proof requires a detailed analysis 
of the crossover regimes from one in which the two layers are effectively decoupled, to 
one where the effective dispersion relation is approximately parabolic around the central 
Fermi points (and the inter-particle interaction is effectively marginal), to the deep in¬ 
frared one, where the effective dispersion relation is approximately conical around each 
Fermi points (and the inter-particle interaction is effectively irrelevant). Such an analy¬ 
sis, in which the influence of the flow of the effective constants in one regime has crucial 
repercussions in lower regimes, is, to our knowledge, novel. 

We expect our proof to be adaptable without substantial efforts to the case where 
74 and A are different from zero, as in (1.5), the intra-layer next-to-nearest neighbor 
hopping 7 q is 0(e), the chemical potential is O(e^), and the temperature is larger than 
(const.)e^. At smaller scales, the Fermi set becomes effectively one-dimensional, which 
thoroughly changes the scaling properties. In particular, the effective inter-particle in¬ 
teraction becomes marginal, again, and its flow tends to grow logarithmically. Pertur¬ 
bative analysis thus breaks down at exponentially small temperatures in e and in U, 
and it should be possible to rigorously control the system down to such temperatures. 
Such an analysis could prove difficult, because it requires fine control on the geometry 
of the Fermi surface, as in [BGM06] and in [FKT04a, FKT04b, FKT04c], where the 
Fermi liquid behavior of a system of interacting electrons was proved, respectively down 
to exponentially small and zero temperatures, under different physical conditions. We 
hope to come back to this issue in the future. 

Another possible extension would be the study of crossover effects on other phys¬ 
ical observables, such as the conductivity, in the spirit of [Mall]. In addition, it would 
be interesting to study the case of three-dimensional Coulomb interactions, which is 
physically interesting in describing clean bilayer graphene samples, i.e. where screening 
effects are supposedly negligible. It may be possible to draw inspiration from the analy¬ 
sis of [GMPIO, GMPllb] to construct the ground state, order by order in renormalized 
perturbation theory. The construction of the theory in the second and third regimes 
would pave the way to understanding the universality of the conductivity in the deep 
infrared, beyond the regime studied in [Mall]. 

Acknowledgments We acknowledge financial support from the ERG Starting Grant 
CoMBoS (grant agreement No. 239694) and the PRIN National Grant Geometric and 
analytic theory of Hamiltonian systems in finite and infinite dimensions. 
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A Computation of the Fermi points 


In this appendix, we prove (3.2). 

-Proposition A.l 


Given 


'V3, 


0 (/c) := 1 + 2e 2 * ^ cos I -^ky 1 , 


the solutions for k G Aqo (see (2.4) and following lines for the definition of A and Aqo) of 


with 


are 




0 < 7173 < 2 


— ^"ikx 


= 0 




Pf,i •= 


1 —7173 
2 


Pf,3 ■— \ 3 


— i arccos 


(A.l) 


P %2 (x i ^^ccos ^ 7173) ^ arccos 

( VI+7173G-7173 ) ^ arccos 


(A.2) 


V 2 JJ 


for CO G {+, -}. 


Proof: We define 


C 




G := 7173 


in terms of which (A.l) becomes 

r 4(2C2 - l)y2 + 2C(2 - G)Y + 1 - G{2C‘^ - 1) = 0 
1 -25(C(4y2 - G) + y(2 - G)) = 0. 


(A.3) 


1 - If 5 = sin((3/2)A:a;) = 0, then kx G {0,27r/3}. Furthermore, since k G Aqo, 
if /ca; = 0 then ky = 0, which is not a solution of (A.l) as long as G < 3. Therefore 
kx = 27r/3, so that G = —1, and Y solves 

4 y2 _ 2(2 _ G)y + 1 - G = 0 
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so that 


Y = 


2-G±G 


and therefore 


, 27r , 2 /1-G 

= ± —^ or = ± arccos - 

^3^/3 ^ ^/3 V 2 


2 - If 5 / 0, then 


G(4Y^ -G) = -Y(2 - G) 


so that the first of (A.3) becomes = 1 + G, which implies 

Y = I ^ n ^ yi + G(2 - G) 

r» 5 I r» 


so that 


, 27r 2 /yrTG(2-G) 

kx = - 1 — arccos - 

3 3 V 2 


or 


, 27r 2 /yrTG(2-G) 

kx = -arccos - 

3 3 \ 2 


, ny 


, ny 


2 

^yi + G 

— 1 ,— cXiL L/GvJo 1 

ys 

1 2 

2 / 

^yiTG' 

— _1_ /— cti 1 

ys \ 

. 2 , 


□ 


B 4x4 matrix inversions 

In this appendix, we give the explicit expression of the determinant and the inverse 
of matrices that have the form of the inverse free propagator. The result is collected in 
the following proposition and corollary, whose proofs are straightforward, brute force, 
computations. 


Given a matrix 


Proposition B.l 


( il a* 0 

a if b 

0 b* ii 

V b 0 c* 


b* \ 
0 


*3 / 


with (j:, 3 ) G and (a, b, c) G C^. We have 


det A = (|b |2 + 3 y )2 + |a|V + \c\Hf + |aP) " 27^e(a*b2c) 


(B.l) 


(B.2) 
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and 


with 




1 

det A 


^ 0a,a 

<6 
S^a 
\ 0^6 


0a,fe 
0a,a 

Qa,b 

0a,a 


0a,a 



0a,a 



0a,& 
0^a 
0a,b 
0a,a y 


Qa,a = -i3|0P - + |cP) 

0a,6=3V-C*((b*)2-a*c) 

0a,a = iia*b + iib*c* 

< 

0a,6 = b((b*)^ - o*c) +3jb* 
0a,6 = -a((0*)^ - a*c) +y2c 

, 0a,a = -*3|aP - n{n + |0P)- 


and given a function 0 (a, b, c,y, 3 ), 

0 +(o, b,c,j:, 3 ) := 0 *(a, b,c,-y,- 3 ). 


(B.3) 


I-Corollary B.2- 

If 3 = y = 0, then 

det^ = jb^ - oc*|% 0. (B.4) 

In particular, A is invertible if and only if b^ 7 ^ ac*. 


C Block diagonalization 

In this appendix, we give the formula for block-diagonalizing 4x4 matrices, which 
is useful to separate the massive block from the massless one. The result is collected in 
the following proposition, whose proof is straightforward. 
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I-Proposition C.l- 

Given a 4 x 4 complex matrix B, which can be written in block-form as 

Bii B^'^ 



(C.l) 


in which B^^, B^'^ and are 2x2 complex matrices and B^^ and are invertible, 
we have 


f 1 0 \ f 1 

V -B^'f>{B^^)-^ 1 y V 0 1 

If B^‘t> - B^'^{B^^)-^B^‘f> is invertible then 


B^^ 0 \ 

0 B^^ - B^^{B^^)-^B^‘^ ) ■ 

(C.2) 


(^00 - B^^{B^^)-^B^‘^)-^ 


is the lower-right block of i? ^ 


D Bound of the propagator in the II-III intermediate regime 

In this appendix, we prove the assertion between (3.38) and (3.39), that is that the 
determinant of the inverse propagator is bounded below by (const.) e® in the intermediate 
regime between the second and third regimes. Using the symmetry under kx e-)• —kx 
and under 27r/3 rotations, we restrict our discussion to cu = -|- and ky —PpQy > 0. 
In a coordinate frame centered at p^q, we denote with some abuse of notation q = 
{kQ,kx,ky) and Pp^ = {0,De^), where e |73 and D = ^^(1 + 0{e^)) (see (3.3)). Note 
that D > 0 IS uniformly bounded away from 0 for e small (recall that 71 = e and 
73 = 0.33e). In these coordinates, we restrict to ky > 0, and the first and third conditions 
in (3.37) read 

^Jkl + e^{kl + kf) ^ -b + {ky - D^Y) ^ ke^, (D.l) 

where kk 2 (J^)^. The second condition in (3.37) implies that {kx + ky) ^ (const.) e^, 
in which case the desired bound (that is, | det ^ (const.) e®, with det A as in (3.38)) 
reads 

81 

e^kl + — \{ikx + ky)'^ - Df{-ikx + ky)f ^ (const.) e®. (D.2) 

Therefore, the desired estimate follows from the following Proposition, which is proved 
below. 
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I-Proposition D.l- 

For all D,e > 0, if (/cq, kx, ky) G satisfies 

ky > 0, \Jkl + e^{kl + > Re^, \Jk^ + e^(9/c2 + [ky - De^Y) > Re^ 

for some constant /t > 0, then, for all a > 0, we have 

f‘kl + a\{ikx + ky)'^ - Df‘{-ikx + ky)^ > Cf', (D.3) 

where 


C := min 



aD"^ a(473 - 3\/I05 )k2 
1 ^’ 288 


T' 


Proof : We rewrite the left side of (D.3) as 

I := e^/cg + a {—k^ + ky — De^ky)^ + ak^ (2ky + . 

If |/co| > then I > K^e®/4 from which (D.3) follows. If |A:o| ^ Ke^/2, then 

kl + kl> 9 / c ^ + {ky - Dff > 

If 1^x1 > (l/4\/3)Ke^, then, using the fact that ky > 0, I > Q;(l/48)D^K^e® from 
which (D.3) follows. If \kx\ ^ (l/4\/3)Ke^, then 


so that 


\ky{ky - De^)\ -kl> ^ R^e^ 


48 


and I > a((3\/l05 — I)^/2304)K^e® from which (D.3) follows. 


□ 


E Symmetries 


In this appendix, we prove that the symmetries listed in (2.32) through (2.38) 
leave ho and V invariant. We first recall 


ho — — 


1 


Xo( 2 -^|A:o|)/I|A| 


5 ( 4" A )()(I) 


/3.L 
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with 


^«(k):=(“« Ai*(k) = A*i{k) 


0 n*{k) 

Q{k) 0 


A'^'^(k) := 


^ ( h>\ /D Sikx 


-f3n{k)e^ 


73^7* (/c)e 


V(V’) = 


(/3|A|)2 




ki-k 2 +k 3 ,a' 


(a,a') ki,k 2 ,k 3 


where 


^,a'{k) ■■= ^ - do,'). 


1 - Global ?7(1). Follows immediately from the fact that there are as many 
as V’~ in ho and V. □ 


2 - 27r/3 rotation. We have 

9.{e^^^^k) = e*'2-^S7(A:), ^ ^-Zih-k 

so that 7j^yl'^'^(T“^k)7k = A'^'^(k) and yl^‘^(T“^k)7k = A^‘^(k). This, together with 
A^^{T~^'k) = A^^(k), implies that ho is invariant under (2.33). 

Furthermore, interpreting as a rotation in around the z axis, 

= da, = d^, e-^"-r-^da = l2 + da, e-*^"2dfe = -Z2 + db, 

which implies, denoting by v{k) the matrix with elements Va,a'{k), 

( Va,a{k) v^i{k) e^^'^^Va,d{h) e-^^'^Ha,b{k) \ 

.^, 3 v-^Jk) viiik) e^^-hv- {k) e-^'^■^Hl^^{k) 

UyC riJ — AU 1 - ^ /,s AU 1- .. /,N . nAU 1 ^ ^ /,N I 


e "-^''-'^VdAk) e "^'’'^v- f^{k) Da,a e ‘^"^'^^Vd,b{k) 

e"’^'^^Vb,aik) e^’^'^'^v^-^ik) e^^’^'^Hb,d{k) Vb,b{h) 


furthermore 


f+ t- 
^ki,a^k 2 ,a 

1 + - 
^ki, 6 ^k 2 ,b 

('^ki'^i )“("^2 '^k 2 )“ 


^ki,a^ki,a 

K,iKd> 
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from which one easily concludes that V is invariant under (2.33). □ 

3 - Complex conjugation. Follows immediately from Q(—k) = and 

v{-k)=v*{k). □ 

4 - Vertical reflection. Follows immediately from 0,{Ryk) = 0,{k) and 

v{Rvk) = v{k) (since the second component of da is 0). □ 


5 - Horizontal reflection. We have ^{Rhk) = ^}*{k), (Tiyl^^(k)iTi = A^^(k), 


cJiA^'^(k)cri = 


0 n{k) 
n*{k) 0 


, cjiyl‘^'^(k)cri = 


iko ^^Q.*{k)e 
73f2(A:)e3*^- iko 


—Sikx 


from which the invariance of ho follows immediately. Furthermore 

^a^a'i.hdhk') ^■7r/i(a),7rh(o') 

where tt/j is the permutation that exchanges a with h and d with b, from which the 
invariance of V follows immediately. □ 


6 - Parity. We have f2(PA:) = n*{k) so that [A^‘^{Pk)]'^ = A«'^(k), [A^^{Pk)]'^ = 
y4'^'^(k), [yl^^(Pk)]^ = yl^^(k). Therefore ho is mapped to 


ho 


Xo{2-^\ko\)m 


^ A««(k) A«^(k) \ 




which is equal to ho since exchanging il) and adds a minus sign. The invariance of V 
follows from the remark that under parity aV’kj a '^P]s .2 a'^^Pki a’ “ ^ 2 ) = 

v{P{k 2 -ki)). □ 


7 - Time inversion. We have 

cj3y4^^(/k)(T3 = — A^^(k), cj3yl^'^(/k)cj3 = —A^^{k), 

a3A't>'f>{Ik)a3 = -A^^{k) 

from which the invariance of ho follows immediately. The invariance of V is trivial. 


□ 
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F Constraints due to the symmetries 

In this appendix we discuss some of the consequences of the symmetries listed in 
section 2.3 on W 2 ^\'k) and its derivatives. 

We recall the definitions of the symmetry transformations from section 2.3: 

• 27 r 

Tk := R^k := {ko,ki,-k2), Rhk := {ko,-ki,k2), 

(h .1) 

Pk := {ko, -ki, -k2), Ik := {-ko, ki, ^ 2 ). 

Furthermore, given a 4 x 4 matrix M whose components are indexed by {a, b, a, b}, we 
denote the sub-matrix with components in {a,b}^ by that with {a,b}^ by 

with {a, b} X {a, b} by and with {a, b} x {a, b} by M'^^. In addition, given a complex 
matrix M, we denote its component-wise complex conjugate by M* (which is not to be 
confused with its adjoint M^). 


I-Proposition F.l- 

Given a 2 x 2 complex matrix M(k) parametrized by k G Boo (we recall that Boo was 
defined above the statement of the Main Theorem in section 1.3) and a pair of points 
{Pf,Pf) e ^lo, if Vk G Boo 

M(k) = M{-k)* = M{R^k) = aiM{Rhk)ai = -a3M{Ik)a3 (F.2) 


and 

Pf = -Pf" = ^-Pf" = RhP% = IPF 
for u) G { —, then C; G '^) £ such that 

^^'(Pf) = dkoM{p%) = fCl, 

M {p%) = va2 , 9fc2 M {p%) = Lowai. 


(F.3) 


(F.4) 


Proof : 

1 - We first prove that M(p‘^) = fiai. We write 

M{p%) =: tl + xai + ya2 + - 20-3 
where {t,x,y,z) G C^. We have 

M{p%) = M(pp“)* = M(p^-) = aiM{p%)ai = -c73M(p^)a3. 
Therefore (f, x, y, z) are independent ofco, t = y = z = 0 and x G M. 
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2 - We now study which we write as 

5fcoM(p‘^) =: tot + xoai + yocr2 + zods. 


We have 

dkoM{p%) = -{dkoM{pp‘^))* = dkoMipp^^) = aidkoM{p%)ai = a3dkoM{p%)a3. 
Therefore {to,xo,yo, zq) are independent of a;, xq = yo = 2:0 = 0 and to G zM. 

3 - We now turn our attention to dk^M: 

dkiM{pp) =: til + xidi + yicj2 + Z1CJ3. 


We have 

dkiM{p%) = -{dk^M{pp‘^))* = dk^M{pp‘^) = -aidkiM{p%)ai = -a3dkiM{p%)a3. 

Therefore {ti,xi,yi, zi) are independent of a;, ti = xi = zi = 0 and yi G M. 

4 - Finally, we consider dkyM: 

dk2M{p%) =: t'f'^1 + x'f ^ai + 2 / 2 "^V 2 + 

We have 

dkiMipp) = -{dk2M{p~‘^))* = -dk^Mipp^) = aidk2M{p‘^)ai = -a3dk2M{p%)a3. 
Therefore t'^'^ = y^'^ = z^^ = 0, G K. □ 

I-Proposition F.2- 

Given a 4 x 4 complex matrix M(k) parametrized by k G Boo and two points (p^, p^) G 
if V(/,/') G {^AY and Vw G {-,+}, 

M//'(p-) = yff'ai, dk,Mff\p%) = 

dk^Mff {p%) = Af'a2, dk2Mff'{p‘^) = uwff'ai 
with , w^^') G independent of tc, and Vk G Boo 

M(k) = M'^(Pk) 

and 

Pf = ^’Pf" 

then 


(F.5) 

(F.6) 

(F.7) 

(F.8) 
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Furthermore, if 
(3/2,-V3/2)) 


then 


= ( 0 ,^, 0 ;^^) and (recalling that 7k = 

M(k) = ( ; ) M(T-k) (; 


e * 72 -fc)a- 3 ^ with I 2 = 
(F.9) 


= -xxj' 






i/?? = = 0 , 

= 0 , = 0 . 


(F.IO) 


Proof : (F.8) is straightforward, so we immediately turn to the proof of (F.IO). 
1 - We first focus on M'^'^ which satisfies 


M'^'^(k) = (F.ll) 

Evaluating this formula at k = p*^, recalling that M'^'^(p‘^) = and noting that 

we obtain = 0. Therefore, deriving (F.ll) with respect to ki, 
i = 1,2, and evaluating at p'^, we get: 

2 

i=i 


with 

r = i(-‘ 

2 I, x/3 -1 

Furthermore, recalling that = z/'^'^cr 2 and 9^2 = ujzu^^ai, 



TpUfciM^^Tp^ = 

which implies 

/ \ 

y ujw'^'^ai j 


1 / coV3{u^^ + \ / CJ2 \ 

4 y -y/3{iJ^‘^ + ro'^'^) io{w^^ - ) \ (Ti ) 


so = —w^ 


2 - We now study which satishes 

M'^«(k) = 


no 



Evaluating this formula and its derivative with respect to ko at h = p'^, we obtain 
Evaluating the derivative of this formula with respect to /c* at k = 

we obtain 

2 

i=i 


Furthermore, 


- ^cj2 +a;^o-i), - a;^cj2), 

which implies 

/ U‘f^^a2 \ 1 / + 3'G7'^^ \ / (72 \ 

y uzu^^ai J 4 y / y <7i y 


so that v'lf’ = 'dj'lf- The case of is completely analogous and gives = 0 

and vf = tuf. 


3 - We finally turn to which satishes 




Therefore for i G {1,2}, 

2 

i=i 

so that = 0, that is = 0. □ 
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